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1 Introduction 

A fairly sophisticated description of electron-phonon superconductivity has existed since the 
early 1960's, following the work of Eliashberg |]l|, Nambu |@], Morel and Anderson 0], and 
Schrieffer et al. |^. All of this work extended the original ideas of Bardeen, Cooper, and 
Schrieffer ||^ on superconductivity, to include dynamical phonon exchange as the root cause 
of the effective attractive interaction between electrons in a metal. For certain supercon- 
ducting materials, Eliashberg theory (as this description is generally called) provides a very 
accurate description of the superconducting state. Nonetheless, as B.T. Matthias was fond 
of iterating 0, this description was never considered (by him and others) particularly help- 
ful for discovering new, high temperature superconductors 0. Part of the problem remains 
that a truly accurate description of the normal state has not been forthcoming. Part of that 
problem is the 'curse' of Fermi Liquid Theory. To the extent that the electron-phonon cou- 
pling causes relatively innocuous corrections to most normal state properties, its underlying 
characteristics remain undetectable (indeed, as will be reviewed here, the characteristics of 
the electron-phonon interaction are made more apparent in the superconducting state). An 
exception may be the A15 compounds, whose anomalous normal state properties might help 
us achieve further understanding of the electron-phonon interaction in these materials 



This review will barely touch upon normal state properties influenced by the electron- 
phonon interaction. A considerable literature continues to develop on this topic, including 
a more microscopic treatment of model systems with simple electon-ion interactions. There 
have been many theoretical developments in the last two decades, many of which have 
been directed towards understanding the high temperature oxides. Some references will be 
provided in the Appendix, but, for the bulk of the chapter, we will focus primarily on the 
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superconducting state in 'conventional' superconductors. In the past, many reviews have 
been written on the role of the electron-phonon interaction in superconductors. The reader 
is directed in particular to the reviews by Carbotte |T^, Rainer |T^, Allen and Mitrovic 



13| , and Scalapino fl^ (they are listed here in inverse chronological order). While we 



have repeated much of what already exists in these reviews, we felt it was important for 
completeness in the present volume, and because the material is presented with a slightly 
different outlook than has been done in the past. 

The first section provides an overview of the subject as we see it, with some details 
relegated to the Appendix. This is followed by a discussion of our knowledge of the electron- 
phonon interaction in metals, including an update on old ideas to use the optical conductivity 
to extract this information. The next two sections provide a very brief review of the impact of 
the electron phonon interaction on the superconducting critical temperature, the energy gap, 
the specific heat, and critical magnetic fields. The next section examines dynamical response 
functions. Again, largely because of the discovery of the high temperature superconductors, 
workers were prompted to re-examine in more detail the effect of stronger electron phonon 
coupling on various response functions. For example, as will be discussed in the pertinent 
subsection, the lack of a coherence peak in the NMR relaxation time was observed. Does this 
(on its own) indicate an exotic mechanism, or can it be explained by damping effects due 
to a substantial electron phonon coupling ? Answers to such questions are reviewed in this 
section. Finally, we end with a summary, including some remarks on various non-cuprate but 
non-conventional superconductors. The Appendix will sketch some derivations and provide 
references to more recent literature. 



2 The Electron-Phonon Interaction: Overview 



2.1 Historical Developments 



The history of superconductivity is an immense and fascinating subject |I5|. While the 
discovery of superconductivity occurred in 1911 [|TB|, from a theoretical point of view, a 
first breakthrough occurred with the discovery of the Meissner-Ochsenfeld effect |]1^, and 
the understanding that this implied that the superconducting state was a thermodynamic 
phase [|l^]. During this time a few attempts were made at proposing a mechanism for 
superconductivity [|^, but, by 1950, when London's book appeared, nothing concerning 



mechanism was really known |2T] . 

In 1950 several important developments took place first, two independent isotope 



effect measurements were performed on Hg [^,|2^, which indicated that the superconduct- 
ing transition was intimately related to the lattice, probably through the electron-phonon 
interaction. These experiments were all the more remarkable because in 1922 Onnes and 
Tuyn had looked for an isotope effect in superconducting Pb, and, within the experimental 
accuracy of the time, had found no effect P3[ . 

Secondly, Frohlich adopted, for the first time, a field-theoretical approach to problems 
in condensed matter. In particular, he studied the electron-phonon interaction in metals, 
and demonstrated, through second order perturbation theory, that electrons exhibit an ef- 
fective attractive interaction through the phonons. Although the theory as formulated was 
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incomplete, it did lay the foundations for subsequent work. In fact one of the essential 
features of this mechanism was summarized in his introduction ||26|: "Nor is it accidental 



that very good conductors do not become superconductors, for the required relatively strong 
interaction between electrons and lattice vibrations gives rise to large normal resistivity." 
His theory correctly produced an isotope effect (recognized in a Note Added in Proof), and, 
moreover, foreshadowed the discovery of the perovskite superconductors, by suggesting that 
the number of free electrons per atom should be reduced. 

After hearing about the isotope effect measurements, Bardeen also formulated a theory 
of superconductivity based on the electron-phonon interaction, wherein he determined the 
ground state energy variationally [^. Both of these theories failed to properly explain su- 



perconductivity, essentially because they focussed on the single-electron self-energies, rather 
than the two-electron instability Another breakthrough occurred a little later when 

Frohlich used a self-consistent method to determine an energy lowering proportional to 
exp(— 1/A), where A is the dimensionless electron-phonon coupling constant. This showed 
how essential singularities could enter the problem, and why no perturbation expansion in 
A would succeed in this problem (although in fact the energy lowering is due to a Peierls 
instability, not superconductivity). 

A parallel development meanwhile had been taking place in the problem of electron 
propagation in polar crystals, i.e. the study of polarons. In fact, this problem dates back 
to at least 1933 when Landau first introduced the idea of a "polarization" cloud due 



to the ions surrounding an electron, which, among other things, renormalized its properties. 



Frohlich also addressed this problem, first in 1937 and then again in 1950 Lee, 
Low and Pines subsequently took up the problem, also using field-theoretic techniques, 
to provide a solution to the intermediate coupling polaron problem. This problem was taken 



on later by Feynman [^], then by Holstein and others [Q, along with many others to the 
present day. In fact, as described in the Appendix, a small group of physicists continues 
to emphasize polaron physics as being critical to high temperature superconductivity in the 
perovskites. 

Pines, having worked with Bohm on electron-electron interactions, and having just used 
field-theoretic techniques in the polaron problem, now combined with Bardeen to derive an 
effective electron-electron interaction, taking into account both electron-electron interactions 
and lattice degrees of freedom [Q. The result was the effective interaction Hamiltonian 



between two electrons with wave vectors k and k' and energies ek and ek' [36|: 



cff 
k,k' 



(k - k')2 + kl 



1 + 



h^uH\i - k') 



;ek-ek')'-^'^'(k-k' 



where is the Thomas- Fermi wave vector, and co'(q) is the dressed phonon frequency. Eq. 
(0) is an effective interaction; a more formal and general approach, utilizing Green functions, 
will be given later. Nonetheless, it is clear that this effective interaction captures the essence 
of "overscreening" , i.e. for electronic energy differences less than the phonon energy, the 
phonon contribution to the screened interaction has the opposite sign from the electronically 
screened interaction, and exceeds it in magnitude. Physically [0, one electron makes a 
transition, which excites a phonon, accompanied by an ionic charge density fluctuation. A 
second electron undergoes a transition caused by this induced charge density fluctuation. If 
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the differences in the electron energies is small compared to the phonon excitation energy, 
the second electron is actually attracted to the first. This is shown pictorially in Fig. 1. 

Eq. d^) represents the starting point for the two-electron interaction in metals. It was fur- 
ther simplified for both the Cooper pair calculation ||3^ and the Bardeen-Cooper-Schrieffer 
(BCS) calculation. The progression of events that ultimately led to a successful theory 
for BCS has been well documented |22|. Most of this part of the story had little to do with 
the details of the attractive mechanism, but rather with the pairing theory itself. Thus, 
one can divide the theory of superconductivity into two separate conquests: first the estab- 
lishment of a pairing formalism, which leads to a superconducting condensate, given some 
attractive particle-particle interaction, and secondly, a mechanism by which two electrons 
might attract one another. BCS, by simplifying the interaction, succeeded in establishing 
the pairing formalism. They were able to explain quite a number of experiments, previ- 
ously performed, in progress at the time of the formulation of the theory, and many that 
were to follow. However, one might well ask to what extent the experiments support the 
electron-phonon mechanism as being responsible for superconductivity [^. Indeed, one of 
the elegant outcomes of the BCS pairing formalism is the universality of various properties; 
at the same time this universality means that the theory really doesn't distinguish one su- 
perconductor from another, and, more seriously, one mechanism from another. Fortunately, 
while many superconductors do display universality, some do not, and these, as it turns out, 
provided very strong support for the electron-phonon mechanism, as initially motivated by 
Frohlich |^ and by Bardeen and Pines |3^. Much of this chapter will be concerned with 
these deviations from universality. 

After the BCS paper appeared, several workers rederived their results using alternative 
formalisms. For example, Anderson used an RPA treatment of the reduced BCS Hamiltonian 
in terms of pseudospin operators and Bogoliubov and others [0^] developed more 
general methods, later to be adapted to inhomogeneous superconductivity by de Gennes 
p3| |. Finally, Gor'kov developed a Green function method, from which both the BCS 
results, and the Ginzburg-Landau phenomenology [|5| could be derived, near the transition 
temperature, Tc. 

The Gor'kov formalism proved to be the most useful, for the purposes of generalizing 
BCS theory (with its model effective interaction) to the case where the electron-phonon 
interaction is properly taken into account in the superconducting state. This was done by 
Eliashberg , as well as Nambu [0] , and later partially by Morel and Anderson and more 
completely by Schrieffer and coworkers [Q^, Around the same time tunneling became 
a very useful spectroscopic probe of the superconducting state besides providing an 
excellent measure of the gap in a superconductor, it also revealed the fine detail of the 
electron-phonon interaction |^9[, to such an extent that tunneling data could be "inverted" 
to tell us about the underlying electron-phonon interactions [0]. These developments have 
been well documented in the Parks treatise |^ . In particular retardation effects are covered 
in the articles by Scalapino [|14| and McMillan and Rowell |5^. An interesting historical 



perspective is provided in the article by Anderson . 

In the meantime, developments in our understanding of the polaron were occurring in 
parallel. The problem of phonon-mediated superconductivity and the problem of the im- 
pact of electron-phonon interactions on a single electron are obviously related, but, after the 
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initial work by Frohlich and Pines and coworkers, the two fields seem to have parted ways. 
Indeed, an excellent summary of the status of polarons at that time is Ref. |Q, where, 
however, there is essentially no "cross-talk" with the theory of superconductivity. Similarly, 
in the treatise by Parks ||5l| there is essentially no discussion of polarons in spite of the 
fact that the 'polaron' really is the essential building block of the BCS theory of supercon- 
ductivity. So, for example, a perusal of the index of the classic texts on superconductivity, 
by Schrieffer [^, Blatt Rickayzen de Gennes |^3[, and Tinkham |^S[ reveals not 
a single entry The reason for this is that the electron-phonon coupling strength in all 
known superconductors was deemed to be sufficiently weak that the only effect on normal 
state properties was a slightly increased electron effective mass. Thus, the electronic state is 
presumed to be well described by Fermi Liquid Theory, upon which the BCS theory (and its 
modifications) is based. It is important to keep this in mind; for this reason we will refrain 
from referring to Eliashberg theory as a strong coupling theory (we ourselves have used this 
term in the past). Eliashberg theory goes beyond BCS theory because it includes retardation 
effects; however, it is still a weak coupling theory, in the sense that the Fermi energy is the 
dominant energy, and the quasiparticle picture remains intact. 

We make this distinction because in recent years polaron theory has experienced a re- 
naissance, and some attempts to explain high temperature superconductivity have utilized 
polaron and bipolaron concepts. The bipolaron is simply a bound state of two polarons, 
analogous to the Cooper pair, except that the latter requires a Fermi sea to exist (at least 
in three dimensions) whereas the former exists as a tightly bound pair in the absence of a 
Fermi sea. In this respect bipolaron theories resemble the quasichemical theory advocated 
by Schafroth and coworkers ||56| , 60|| in the 1950's. Tightly bound electron pairs are now rec- 
ognized as the strong coupling limit of the BCS ground state; the transition to the normal 
state is, however, governed by very different (and as yet undetermined) excitations compared 
to BCS theory. We will refer to some of this work in the course of this chapter. 

To complete this brief historical tour, we should add that in 1964, with the suggestion of 
a theorist |^, what has emerged as a new class of superconductors was discovered [Q. The 
actual superconducting compound was doped Strontium Titanate (SrTiOs), a perovskite 
with low carrier density. This compound, along with BaPbo.75Bio.25O3, another doped per- 



ovskite discovered in 1975 ||63| with a transition temperature of 12 K, were the precursors 



to the modern high temperature superconductors discovered by Bednorz and Miiller 



In fact, with fortuitous foresight, Schooley et al. ||6^ remarked, "If SrTiOa had magnetic 
properties, a complete study of this material would require a thorough knowledge of all of 
solid state physics." Little did they know that in 1986 perovskites would be discovered, that 
not only had high superconducting transition temperatures, but also exhibited a plethora 
of magnetic phenomena. We should also note that the so-called cuprates, which presently 
exhibit superconducting transition temperatures up to 160 K (under pressure), all contain 
Cu02 layers, whereas the cubic oxides (such as SrTiOs, BaPbo.75Bio.25O3, and Bai_x.K^.Bi03 



1^ (with ^ 30 K)) do not. For this reason many workers have come to regard the layered 



cuprates and the cubic oxides as belonging to two completely separate (and unconventional) 
classes, even though they are both essentially low carrier density perovskites. 
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2.2 Electron-Ion Interaction 



2.2.1 Overview 

A useful ab initio theory has to begin from some fundamental starting point. In condensed 
matter systems the starting point is usually taken to be electrons and ions (with their 
charges, and masses, etc.) along with the chemical composition of the material ||12|. Given 
these ingredients, the prescription for calculation is, in principle, straightforward. One has 
to solve the many-body Schrodinger equation, with a Hamiltonian consisting of one-body 
kinetic energy terms and the two-body Coulomb interaction. The form of these terms, along 
with all the constants involved, are known, so all that is required to solve the problem is 
perhaps some ingenuity along with unlimited computer resources. This has been referred to 
by Laughlin as the Condensed Matter version of "The Theory of Everything" . 



Of course the difficulty is that, even if one could solve this problem, one would not 
recognize what the solution represented. The notion of ionic collective modes (i.e. phonons), 
for example, would not be very transparent in such an approach. More obscure still would 
be the distinction between a superconducting state versus a metallic state. 

Instead, an approach which separates the complex many-body problem into smaller, more 
tractable pieces, has traditionally been adopted in condensed matter, and in particular in the 
problem of superconductivity p,|12|,|14|. The most systematic approach has been discussed 
by Rainer [|12| . The premise in this approach is the observation that many metals (amongst 



which many undergo a transition to a superconducting state) are well described by Landau 
Fermi Liquid Theory. This allows for an asymptotic expansion in small parameters like 
ksTc/Ep, hujphon/ Ep and 1/kpi, where Ep {kp) is the Fermi energy (wavevector) , cjphon is a 
typical phonon frequency, and £ is the electron mean free path. He separates the problem into 
the "/izg'/i energy problem" (effect of Coulomb interactions amongst the electrons themselves 
as well as between the electrons and the fixed nuclear potentials), and the "/ow energy 
problem" (the dressing of conduction electrons with phonons), and the eventual formation 
of the superconducting state. Most of this review will concern the low energy problem. 
In our opinion the high energy problem is not at all solved at present, from a truly "ab 
initio" approach. For example, strictly speaking, one cannot rely on any of the expansion 
parameters mentioned above, because one does not know, in principle, whether one has 
a metal with a well-defined Fermi surface, to begin with. Nonetheless, by appealing to 
experimental observation, one can use for many cases the fact that nature has already solved 
the high energy problem, and proceed from there to solve the low energy part. This has been 
the dominant philosophy throughout most of the last four decades towards understanding 
superconductivity. 

The difficulty with this approach was exemplified by the discovery of superconductivity 
in the layered perovskites; band structure calculations for the parent compound (La2Cu04) 
demonstrated that it was a metal, when in fact the real material was an antiferromagnetic 
insulator. This problem was later repaired ^7j, but it remains the case that band structure 



calculations fail to properly take into account strong Coulomb correlations, and remain 
somewhat powerless to reliably predict a breakdown of the Fermi Liquid picture. 

With these caveats, the "ab initio" approach of Ref. [0] has experienced excellent success 
in cases where a metallic state is known to exist, and experimental input has been used in 
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the theory. We will comment in particular on the "low energy" part of the theory later in 
this chapter. A thorough discussion is available in Ref. . 



2.2.2 Models 

The net result of a proper handling of the "high energy" problem in the case of a well- 
behaved metal is a set of input parameters for the low energy problem that are simple enough 
to make the remaining part of the problem appear to have arisen from a non-interacting 
model. The distinction is that the input parameters (band structure, phonon spectrum, 
etc.) come not directly from specified model parameters, but rather from previous calculation 
and/or experiment. For this reason, we now discuss possible models for the electron-phonon 
interaction, which, for the moment, we view as fundamental models in their own right, and 
not as models which somehow parameterize (and disguise) the "high energy problem". 

The reason for this is that we hope to accomplish several tasks simultaneously. First, we 
will in effect work through the "low energy problem" discussed in the previous subsection. 
Secondly, we will touch upon some of the more recent work on electron-phonon Hamiltonians, 
which are characterized not so much by comparison with experiment as comparison with 
some "exact" solution, as attained, for example, by Quantum Monte Carlo methods |]68|,|69|. 
Thirdly, we will also be able to make contact with recent ongoing work on the polaron (and 
bipolaron). These latter two topics are presented here more by way of a digression. Some 
further detail is presented in an Appendix, but for a more thorough discussion the cited 
literature will have to be consulted. 

It is always tempting to immediately compare the results of a calculation with experiment; 
agreement justifies the starting model (in this context this would mean the Hamiltonian, with 
associated parameters), whereas disagreement would tend to rule out the starting model as 
a candidate. In the many-body problem, however, life is not so simple. For one thing, 
we know the starting Hamiltonian, as emphasized in the previous subsection. We will get 
agreement with experiment if we were only able to routinely calculate any observable. How- 
ever, in our endeavour to understand many-body systems, we have grown to utilize effective 
Hamiltonians, which would capture the essence of the phenomenon under investigation. The 
purpose of this strategy is twofold; we make sense of the many-body system in terms we can 
understand, and we make the calculation itself more tractable in practice. 

There are many Hamiltonians in condensed matter physics, which were derived as effec- 
tive Hamiltonians for some particular problem, but, which have since taken on a life of their 
own. This is true because (a) they have withstood solution in spite of their simplicity, and 
(b) they epitomize some qualitative aspect of the more general problem. Famous examples 
are the Heisenberg/Ising model for spins, and the Hubbard model for fermions with spin 
degrees of freedom. In the electron-phonon problem several effective models have arisen over 
the years, the three most prominent of which have been the Frohlich Hamiltonian the 
Holstein model [Q, and the BLF (Barisic-Labbe-Friedel) model |]70[ (also known as the SSH 



(Su-Schrieffer-Heeger) model [0). The Frohlich Hamiltonian was derived in a continuum 
approximation (see Ref. [|75| or [|7^] for a derivation), and results in a coupling between 
the electron density and the ionic momentum (a canonical transformation changes this to 
the ionic displacement) which diverges as the momentum transfer between electron and ions 
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goes to zero. This Hamiltonian has been the subject of many investigations of the polaron. 
Holstein proposed his model as a simphfication in which the interaction between electron and 
ion is more local; in fact in some ways the simplification Hubbard WM invoked to replace the 



long-range Coulomb interaction is analogous to the simplification that the Holstein model 
represents compared to the Frohlich Hamiltonian. Both the Frohlich and Holstein models 
represent couplings of the electron to an optical phonon mode. We will focus on the Hol- 
stein model since it is particularly amenable to numerical simulations. In contrast, the BLF 
(SSH) model couples the electron to the relative displacement of nearby ions, i.e. an acoustic 
phonon mode. The physics is simple; in the Holstein model ionic distortions affect the elec- 
tron energy level at a particular site, while in the BLF model ionic displacements affect the 
electron hopping amplitude. These are represented pictorially in Fig. 2, although of course 
the coupling is dynamic. 

The BLF model gained prominence in the 1980's when it was used to describe 
solitons in conducting polymers; otherwise comparatively little effort has been expended 
towards an understanding of its properties, particularly in two or three dimensions. The 
BLF Hamiltonian is 

^ = E — + H -Kin, - u,f 

- E ~ " ■ ("i ~ "i)(4<TCia + h.c), (2) 

<ij> 
a 

where the first line refers to the ions, with mass M and spring constant K. The ionic 
degrees of freedom are described by the ion momentum, pj, and displacement, Uj, at site i. 
The electrons are described by creation (annihilation) operators cl^. (qo-) for an electron with 
spin (J at site i. The electron hopping amplitude is given by tj-,-; this in turn is modulated 
by ionic vibrations, and therefore results in the electron- ion coupling with strength The 
coupling constant |a| is proportional to the gradient of the hopping overlap integral between 
electron orbitals on two neighbouring sites. 

Equation ^ gives rise to the standard electron-phonon Hamiltonian, as written in mo- 
mentum space: 

^ = EekcLck<x + ^/i^qa^jflq + -^5I^(k>k')K-k' + a^_(k_k/))cj,vCka. (3) 



kcr q V kk' 

(7 

We have used the conventional oscillator operators, Oq = ^^^^(uq + ipq) and the standard 
Fourier expansions, cl^. = X^k c*'*''^'cJ.o-, etc. The phonon dispersion is given by cjq, where, 
in principle, q includes branch indices as well as momenta within the first Brillouin zone, 
and 5f(k, k') is the coupling function. For the BLF Hamiltonian, this coupling function has 
a very specific form (involving sine functions). A more general consideration of the electron- 



ion interaction yields a Hamiltonian of essentially the same form [0, |T^, but where the 
parameters involved are understood to already contain the "high energy" effects alluded 
to earlier. State-of-the-art computations of the electron-ion coupling strength, are given, 
for example, in Ref. (for La2-xSr2,Cu04) and in Ref. (and references therein, for 
AsCeo)- 
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The Holstein Hamiltonian is 

H=-tY. {clcj, + H.c.) + Y^i^h + ^Kxj] -aY^x^n,^, (4) 

<i3> i ^^^^ ^ ia 

a 

where the parameters are as before except that the displacement variable Xi represents the 
(one-dimensional) displacement of some optical mode (say a breathing mode) associated 
with the iih site, and the electron-ion coupling a represents the change in site energy (per 
unit displacement) associated with this mode. In momentum space this Hamiltonian is 
particularly simple: 

i^ = I]ekcLcka + Y.^^Ea\a^ + ^ + «-(q))4+qaCka, (5) 

k<7 q ViV kq 



where 0;^; is the Einstein mode frequency and g = y " . This model has been studied 
extensively in the last twenty years, at least partly due to its simplicity. Some of this work 
is reviewed in the Appendix. 



2.3 Migdal Theory 

The primary language of many-body systems is the Green function, or propagator. Many 
books have been written (see for example Refs. [^8^) about the Green function formalism, 
so we will bypass a thorough discussion here. A sketch of the derivation of the Migdal |^ 
equation for the electron self-energy is given in the Appendix. Migdal argued that all vertex 
corrections are 0(m/M)^/^ compared to the bare vertex, and therefore can be ignored. Here 
m (M) is the electron (ion) mass. This represents a tremendous simplification, and allows 
one to solve a theory which should work for arbitrary coupling strength (this is, in fact, not 
the case, for reasons that will become apparent in the next section). 

An "exact" formulation of the electron-phonon problem can be summarized |]84|-p^ in 
terms of the Dyson equations (written in momentum and imaginary frequency space): 

G(k, iujm) = [Go(k, iujm)~^ - S(k, ium)]'^ (6) 

for the electron, and 

D{q, iVn) = [^o(q, ii^ny^ - n(q, zz/„)]"^ (7) 

for the phonon, where G(k, iuJm) is the one-electron Green function, /^(q, «//„) is the phonon 
propagator, and E(k, itUm) is the electron and H(q, it/„) the phonon self energy. Then, 

E(k, iuJm) = -^-3 5'k,k'-D(k - k', iuJm-iuJm')G(k' , iuJm')T(}i, iu^rn'] k, iuJm] k - k', iuJrn-iuJn 

^ k',m' 

(8) 

and 

2 

H(q, iun) = — Y 5'k,k+qG'(k + q, ium + iUn)G(k, iuJm)T(k + q, ium + i^n] k, iUm] q, iJ^n), 

k,m 



(9) 



10 



where the vertex function F can only be defined in terms of an infinite set of diagrams (i.e. 
not in closed form). 

The non-interacting propagators are 

Go(k, iujm) = [i(^m - (ck - (10) 

for the electron and 

D,iq,tu^) = [-M{uj\q) + i^l)]-' (11) 

for the phonon, where ek is the single electron dispersion (band indices are implicit here and 
in the following), fi is the chemical potential, and uj{q) is the phonon dispersion. In writing 
these relations we have adopted the finite temperature Matsubara formalism, with Fermion 
{iUm = mT{2m — 1)) and Boson (iz/„ = i2T:Tn) Matsubara frequencies, where m and n are 
integers and T is the temperature {ks = 1). The Matsubara sums in Eqs. (^,^ extend over 
all integers, and the momentum sums extend over the first Brillouin zone. This convention 
will be maintained unless noted otherwise. 

Migdal's approximation was to set the vertex function F equal to the bare vertex, g. 
Then, the electron self-energy can be written: 

j:{k,iujm) = |fi'k,k'P^(k - k',ia;m - iu;m/)G(k',iu;m/). (12) 



Migdal [0] also included renormalization effects in the phonon propagator. With an appli- 
cation to real materials in mind, however, the electron dispersion relations will have been 
obtained from a band structure calculation, and the phonon properties will generally have 
been taken from experiment. In this case the phonon self energy is omitted entirely (to avoid 
double counting). In addition electron-electron effects have been omitted, as they have been 
presumed to be included already in the band structure and phonon calculations (to the best 
extent possible). 

Alternatively, Eq. (|12[) can be viewed as having been derived from some microscopic 
electron- ion Hamiltonian. For example, in the case of the Holstein Hamiltonian, Eq. (^), 
S'k.k' — ^ 9, the constant appearing in Eq. (|), and the electron band structure is given by 
Ck = —2tcos{kx) (in one dimension, and for nearest-neighbour hopping only). In addition, 
the phonon frequency becomes dispersionless (co'(q) ue) and the phonon self energy is 
given by some appropriate approximation. Such an identification is useful for comparison to 
exact results (usually done numerically - see the Appendix for references). 

In the classical hterature [|^,|8^,|8^,|8^,|88|, Eq. (|12D is simplified in the following way. 
First, very often the phonon propagator is provided separately, usually by inelastic neutron 
scattering measurements [39,^|. To see how, one first writes the phonon propagator in 



terms of its spectral representation [jT3|: 



2u 

D{q,iuri) = duB{q,u)- — (13) 



[tiy„ y - z/2 



where S(q, v) is the phonon spectral function 



S(q,z/) = -ilmD(q,z/ + Z(5). (14) 

TT 
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The spectral function is positive definite, and obeys a sum rule; it is the quantity that 
is constructed with fits to high-symmetry phonon dispersion curves measured by inelastic 



neutron scattering Following this tact a calculation of the phonon self energy is no 

longer required. Another simplification was recognized in Ref. |^5|; this is the use of the 
non- interacting electron Green function Go(k, iujm) in the right hand side of Eq. (|I2|) instead 
of the full self-consistent choice, G(^,iuJm)- This approximation is valid when particle- 
hole symmetry is present and the infinite bandwidth approximation is invoked. This latter 
approximation is used extensively in the early literature on metals and superconductors; a 



systematic explanation of the logic is provided in Ref. ITH], and requires the usual hierarchy 
of energy scales, cUphon << Ep {h = 1). The result is 

S(k,zu;^) = — ^ / rfi/|5(k,k'r5(k-k',z/)- — — -Goi}^ .iuj^')- (15) 

The form of Eq. ([T5|) allows one to introduce the electron-phonon spectral function, 

a2F(k,k',z/) =iV(/i)|(7k,kf 5(k-k',z/), (16) 
where N{^) is the electron density of states at the chemical potential. At this point one 



can introduce 'Fermi surface Harmonics' [13, 91], and define an electron self-energy with 



Fermi momentum which depends on Matsubara frequency, and on the angle around the 
Fermi surface. Elastic impurities would act to homogenize the self-energy (as well as other 
properties), so a more useful function for dirty superconductors is the Fermi-surface-averaged 
spectral function, 

a'F{v) ^ -i- «'^(k, k', v)5{e^ - f^)6{e^, - /x). (17) 

yt^) k,k' 

To gain an understanding of electron-phonon effects, Englesberg and Schrieffer solved 
this model for two simple phonon models, the Einstein and Debye models. Here we summa- 
rize their results for the Einstein model, with unmodified phonon spectrum, a simpler case 
since both the phonon spectrum and the bare vertex function are independent of momentum. 
In this case ^^k.k' = 9 and -B(q, z/) = 5(i/ — uje)- Using, in addition, the prescription 



k 



along with a constant density of states approximation, extended over an infinite bandwidth, 
one obtains for the electron self energy 

/°° 1 . 1 1 
de-^Y. 2 I / ? V' (19) 
-oo P ^ Uf^ + [UJrn' - ^m)'' i(^m' - - N 

where we have used the standard definition for the electron-phonon mass enhancement pa- 
rameter, A: 

X^^Td.'^, (20) 

JO y 
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which, for the Einstein spectrum used here, reduces to 

A = 2N{eF)gyuE. 
Performing the Matsubara sum yields 
Xuje 



de 



uje - (e- yU) iujm + CJ^; - (e - /i) 



(21) 



(22) 



where /(e — fi) is the Fermi function and n{ujE) is the Bose distribution function. The 
remaining integral can also be performed 



S(z) 



Xujf 



-2miniLUE) + 1/2) + ^(- + i'^^^^) - ^(- - t'^^^^] 



(23) 



where il}[x) is the digamma function [|13|,^ and the entire expression has been analytically 
continued to a general complex frequency z. Because we performed the Matsubara sum first, 
before replacing iiOm with 2;, this is the physically correct analytic continuation . 



At zero temperature one can use well-documented properties of the digamma function, 
or, more simply, refer to the analytic continuation of Eq. (122|), since the Bose and Fermi 
functions may be more famihar. Since u^uje) — and /(e — fi) ^ 9{fi — e) as T — > {9{x) 
is the Heaviside step function), the self energy at T = is 



XuJr 



In 



Ue- z 



\UJe + z 



(24) 



Spectroscopic measurements yield properties as a function of real frequency; because of the 
analytic properties of the Green function, this corresponds to a frequency either slightly 
above or below the real axis. We will use frequencies slightly above, and designate the 
infinitesmal positive imaginary part by Thus, 



E(cj + 16) 



XuJi 



In 



Ue 



UJe + 00 



-mt 



-UJe 



(25) 



The real and imaginary parts of this self energy are shown in Fig. 3, along with the non- 
interacting inverse Green function {u — (ck — /i)) to determine the poles of the electron Green 
function (see Eq. (|])) graphically. A quantity often measured in single particle spectroscopies 
is the spectral function, A(k, u) defined by 



A{k, u) = --^^(k, uj + i6). 
With this definition, we obtain, through Eq. 

A{k,uj) = 6{uj 



(ek - /i) - ^ In , 

Xuje/2 



) and 
UJe ~ ^ 



UJe + UJ 



UJ 



(ek-/i)-^ln 









1 uje~^ I 




/ ttXuje \ 


1 loe-\-^ ' 




\ 2 ; 



(26) 

if I |< UJe, 

if I |> ue- (27) 
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Plots are shown in Fig. 4. Each spectral function displays a quasiparticle peak, whose 
strength Ok and frequency cuk is implicitly dependent on wavevector 



ak = (1+ . (28) 

where tUk is the solution (between —oje and oje) to the zero of the delta-function argument 
in Eq. (|27D . For all momenta (or equivalently all Ck ~/^) there is a solution, whose frequency 
approaches uje asymtotically as ek — /i ^ 00. The weight of this peak starts at the Fermi 
surface (ek = /i) as 1 /(1+A) and quickly goes to zero according to Eq. (p8|) as tUk uje, which 
occurs for ek'^2ujE- For larger ek a quasiparticle peak forms once again, albeit with non- 
zero width, at approximately the non-interacting electron energy, ek = H- At intermediate 
ek ~ uje, the quasiparticle picture has broken down, and a description as described here is 
required for a complete picture. 

How well the Migdal approximation works in specific circumstances is the subject of on- 
going research (see, for example, Refs. p^-p8|, and the Appendix. For example, Alexandrov 



et al. found an apparent breakdown (for coupling strengths greater than 1, within the 
Holstein model) to the approximation when a finite electronic bandwidth was taken into 
account. 

We have focussed on the modifications to the electron spectral function due to the 
electron-phonon interaction. For excitations at the Fermi level (ek = /i), the quasiparti- 
cle pole remains there (i^kp = 0)^ remains infinitely long-lived (it is a delta-function), but 
has a reduced weight, by a factor of 1 + A. This same factor enhances the effective mass. 



and alters various normal state properties in a similar way [p8| , |100|| . For example, the low 
temperature electronic specific heat is linear in temperature with coefficient usually denoted 
by 7, which is proportional to the electron density of states. The electron-phonon interaction 
enhances this coefficient by the same factor, 1 + A. Other renormalizations are reviewed in 



Ref. 88 



2.4 Eliashberg Theory 

Eliashberg theory is the natural development of BCS theory to include retardation effects due 
to the 'sluggishness' of the phonon response. In fact, insofar as BCS introduced an energy 
cutoff, ujd (the Debye frequency), they included, in the most minimal way, retardation effects. 
However, Eliashberg theory goes well beyond this approximation, and handles momentum 
cutoffs and frequency cutoffs separately. We begin this section with a very brief review of 
BCS theory, followed by a more detailed discussion of Eliashberg theory. 



2.4.1 BCS Theory 

Before one establishes a theory of superconductivity, one requires a satisfactory theory of 
the normal state. In conventional superconductors, Fermi Liquid Theory appears to work 
very well, so that, while we cannot solve the problem of electrons interacting through the 
Coulomb interaction, experiment tells us that Coulomb interactions give rise to well-defined 
quasiparticles, i.e. a set of excitations which are in one-to-one correspondence with those 
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of the free-electron gas. The net result is that one begins the problem with a 'reduced' 
Hamiltonian, 

-f^red = XI ^kcLckfT + ^k,k'Ck|clkjC_k'iCkt , (29) 

kff kk' 

where, for example, the electron energy dispersion already contains much of the effect due 
to Coulomb interactions. The important point is that well-defined quasiparticles with a well- 
defined energy dispersion near the Fermi surface are assumed to exist, and are summarized 
by the dispersion e^. The pairing interaction VCk, k') is assumed to be 'left-over' from the 
main part of the Coulomb interaction, and this is the part that BCS simply modelled, based 
on earlier work by Frohlich |2^ and Bardeen and Pines [p5|] . 

Complete derivations of BCS theory have been provided elsewhere in this volume; here 



we state the final result 



-E^kk'^tanh^, (30) 



where 



E^ = Jie^- f^y + Ai (31) 



is the quasiparticle energy in the superconducting state, and Ak is the variational parameter 
used by BCS. An additional equation which must be considered alongside the gap equation 
(0) is the number equation, 

n=l-lv^^^tanh^. (32) 

Given a pair potential and an electron density, one has to 'invert' these equations to determine 
the variational parameter Ak and the chemical potential. For conventional superconductors 
the chemical potential hardly changes on going from the normal to the superconducting 
state, and the variational parameter is much smaller than the chemical potential, with the 
result that the second equation was usually ignored. 

BCS then modelled the pairing interaction as a negative (and therefore attractive) con- 
stant with a sharp cutoff in momentum space: 

V^k,k' ~ -Ve{ujD- I (ek - f^) mujn- \ (ek' - /x) I). (33) 

Using this potential in Eq. (|30D , along with a constant density of states assumption over the 
entire range of integration, we obtain 

1 f^D de BE , 

where A = N{fi)V. At T = 0, the integral can be done analytically to give 

A = . (35) 

1 — exp (— 1/A) 
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In weak coupling this becomes the more famihar 



A = 2cJz5exp(-l/A), (36) 

while in strong coupling we obtain 

A = 2ujdX. (37) 



Both of these results are within the realm of BCS theory (at zero temperature) il0l| , [l02| 



although the latter generally requires a self-consistent solution with the number equation, 
Eq. (H). 

Close to the critical temperature, T^, the BCS equation becomes 

1 /•/3'^d/2 tanhx 

ax , (38) 



\ Jo X 

which can't be solved in terms of elementary functions for arbitrary coupling strength. 
Nonetheless, in weak coupling, one obtains 

Te = 1.13cJi5exp(-l/A), (39) 

and in strong coupling 

T, = ujdX/2. (40) 

It is clear that or the zero temperature variational parameter A depend on material prop- 
erties such as the phonon spectrum (ujd), the electronic structure (iV(/i)) and the electron-ion 
coupling strength {V). However, it is possible to form various thermodynamic ratios, which 
turn out to be independent of material parameters. The obvious example from the preceding 
equations is the ratio In weak coupling (most relevant for conventional superconduc- 

tors), for example, we obtain 

2A 

Tl^ = 3-53, (41) 

a universal result, independent of the material involved. Many other such ratios can be 
determined within BCS theory, and the observed deviations from these universal values 
contributed to the need for an improved formulation of BCS theory. For example, the 
observed value of this ratio in superconducting Pb was closer to 4.5, a result that is readily 
understood with Eliashberg theory. It is worth noting that simply extending BCS theory to 
the strong coupling limit (see Eqs. (|37| , |40|) above) results again in a universal constant, — 



4, which is the maximum value attainable within BCS theory with a constant interaction 
103|| , and is still clearly too low. 



Other aspects of BCS theory, particularly those which prove to inadequately account for 
the superconducting properties of some materials (notably Pb and Hg) will not be reviewed 
here. Instead, we will make reference to the BCS limit as we encounter various properties 
within the experimental or Eliashberg context. 
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2.4.2 Eliashberg Equations 



In most reviews and texts tliat derive the Eliashberg equations, the starting point is the 
Nambu formahsm ||^. While this formalism simplifies the actual derivation, it also provides 
a roadblock to further understanding for the uninitiated. For this reason we have followed the 
conceptually much more straightforward approach (provided by Rickayzen [^, for example) 
in the derivation outlined in the Appendix. The result can be summarized by the following 
set of equations: 



S(k, ium) 



0(k, iujj; 

G(k, iun 
F(k, iur, 



— y 



F(k',2t^^/), 



^-^(k, it^^)G-i(-k, -iujyn) + 0(k, iu;m)(^(k, ia;^) 

0(k, iujm) 

^-^(k, 2u;^)G'-i(-k, -iuJm) + 0(-k, -^a;^)0(-k, -iUm) 
G<7^(k,iu;„) - S(k,it^„). 



(42) 
(43) 

(44) 

(45) 
(46) 



Another couple of equations identical to Eqs. (|^) and (145|) , except with and F instead of 
and F, have been omitted; they indicate that some choice of phase is possible, which will 
be important for Josephson effects ||104|| but not for what will be considered in the remainder 
of this chapter. Therefore, we use = ||105|| . 

Note that (^"^(k, -icj^) is the inverse of the non-interacting Green function, in which 
Hartree-Fock contributions from both the electron-ion and electron-electron interactions are 
assumed to be contained. 

Following the standard practice we have used a kernel given by 



Akk'(^) 



dv 



(47) 



where a^T^Fiv) is given by Eq. (|TB|). Eqs. ( P^ - ^T] ) have been written in a fairly general 
way; in this way they can be viewed as having arisen from a microscopic Hamiltonian as in 
Eqs. (0-^) (although electron-electron interactions have been included in the pairing channel 
only, and not in the single electron self energy), or, alternatively, from a treatment of real 
metals, where, as mentioned earlier, the electron and phonon structure come from previous 
calculations and/or experiments. These equations emphasize the electron- ion interaction; 
attempts to explain superconductivity through the electron-electron interactions have been 
proposed in the past, mainly through collective modes [ 106| - 113 |; some of these attempts will 
be treated elsewhere in this volume in the context of high temperature superconductivity. 

Assuming the electron and phonon structure is given, Eqs. ( ^^ - |47| ) must be solved for the 
two functions, S(k, icUm) and 0(k, icUm)- The procedure is as follows: it is standard practice 
to separate the self energy, E(k, -idj^), into its even and odd components W^ '- 



[S(k, ium) - S(k, -iu^)] 
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xiKi^^m) = ^[T.{k,iuj^) + T,{k,-iuj^)] (48) 

where Z and x ^i-re both even functions of iu^ (and, as we've assumed all along, k). Then, 
Eq. (^21) becomes two equations. 



Z(k ' ) - 1 + — V -^'^'"'(^^^ "~ ^^'"'^ {uJm'/l^m)Z{}i',iuJrn') .^gx 

Npj^, iV(/i) u;^,Z2(k',za;„O + (ek'-/i + x(k',«cu„O)' + 02(k',^u;„O 

^/^kw ^('") u;^,Z2(k',za;„O + (ek'-/^ + x(k',^^™O)' + 0'(k',^^„.') 
along with the gap equation (Eq. (^)): 

0(k,^a;^) = — ^ -— Vkk' 



k^' ^ ^(/^) ' ''''' ^ ^^'^'(k', tuo^') + (ek' - /i + X(k', iuj^^)f + 02(k', ^o;, 

(51) 

These are supplemented with the electron number equation, which determines the chemical 
potential, /x: 



^ = 1^ E G(k',zu;^Oe^"'"'°'' (52) 



2 

^kw 

%^m') + (ek' - + x((k', ^a;^0)' + 0'(k', ^^n^') ' 

These constitute general Eliashberg equations for the electron-phonon interaction, in which 
electron-electron interactions enter explicitly only in the pairing equation. Very complete 
calculations of these functions (linearized, for the calculation of T^) were carried out for Nb 
by Peter et al. [|114|| , and for Pb by Daams ||115|| . 

The more standard practice is to essentially confine all electronic properties to the Fermi 
surface; then only the anisotropy of the various functions need be considered. Often these 
are simply averaged over (due to impurities, for example), or the anisotropy may be very 
weak and therefore neglected. In this case the equations (|49| - |53D can be written 



Zm = 1 + 7rrVA(iu;m -2t^m')-7=^==^^=^o("^') (54) 



Xm = -7rT^X{iuj„i- iuJrn')Ai{m') (55) 

m' 

(j)m = TrT^(x{iuJrn- i(^m') - N{n)Vcoui] I ^^"^ ^ ^o("^0 (56) 

l-27rTA^(/i)^Ai(m') (57) 



n 



where we have adopted the shorthand Z{iujm) = Zm-, etc, \{z) and Koui represent appropriate 
Fermi surface averages of the quantities involved, and the functions y4o(m') and y4i(m') are 
given by integrals over appropriate density of states, using the prescription (0) to convert 
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from Eqs. (^9| - |53[) to Eqs. (|5^^7D . If the electron density of states is assumed to be constant, 
then, with the additional approximation of infinite bandwidth, ylo(m') = 1 (actually a cutoff, 
0{uJc— I uJm' I), is required in Eq. (|56|)), and Ai{m') = 0. This last result effectively removes 



Xm (and Eqs. ([55|j57|) ) from further consideration. An earlier review by one of us ||TT 
covered the consequences of the remaining two coupled equations in great detail. 

Nonetheless, a considerable effort has been devoted to examining gap anisotropy, as well 
as variations in the electronic density of states near the Fermi surface. We describe some of 
this work in the following few paragraphs. 

Referring back to Eqs. ( ^5| - |5BD , one can rewrite the summation over k' on the right- 
hand-side of these equations as an integral over energy plus an integral over angle (for a 
given constant energy surface). In carrying out the energy integration the energy dependent 
electron density of states (EDOS), A^(e), introduces a new weighting factor if iV(e) exhibits 
variations over the energy scale of the phonon frequencies. On the other hand, the integration 
over angle will account for variations of the gap and other quantities in the integrands 
with momentum direction. There is a large literature on each of these complicating effects, 
starting with anisotropy effects | [Ll(j| , |117| , and more recently with EDOS energy dependence 



T3|,[TT8Hl20f 



Concerning anisotropy, the observed universal decrease in T^. with increasing impurity 
concentration (i.e. so-called 'normal' impurities, deemed to be innocuous by Anderson's 
arg ument |12T| ) can be attributed to the washing out of gap anisotropy. To see why this 
decreases Tc (we omit here effects due to valence changes) we note that the impurity potential 
scattering has a tendency to homogenize the gap on the Fermi surface. This tends to reduce 
the gap in some directions, and it is these directions that make the maximum contribution 
to Tc, and so is reduced. A simple BCS calculation can demonstrate this analytically. 
One makes a separable approximation for the pairing potential, Eq. (p3|), to be used in the 
BCS equation (^): 

^k,k' = -^(l + ak)(l + ak'), (58) 

where the same energy cutoffs are assumed, and is a function of momentum direction 
only. Assuming Ok to be small with a Fermi surface average equal to zero (i.e. < Ok >= 0) 
and Ok = a^, with <> denoting an angular average over the Fermi surface, then clearly 
Ak = Ao(l + flk)- Solving the resulting equation yields 

< Ak >= Ao = 2a;^exp(- ^^^|^^^ )(l - ^a^) (59) 

in the weak coupling approximation. Similarly, one can solve the Tc equation, to obtain 

Tc = l.lSo^D exp {- ] ). (60) 

This last equation demonstrates that Tc is increased by anisotropy. Hence, increased scat- 
tering due to impurities will decrease Tc, as the anisotropy is washed out. Finally, the gap 
ratio, 

^^3.53(.-|.^), («) 
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showing that anisotropy reduces this quantity. 

How big can the anisotropy be in pure conventional superconductors ? Microscopically 
the anisotropy is related to band structure anisotropy plus anisotropy in the electron-phonon 
spectral function from Eq. (p!6|) , a^F(k, k', u). In Fig. 5 we show the results of a calculation 
of the gap anisotropy in Pb as a function of position on the Fermi surface |p.22||. These 



calculations include multiple-plane-wave effects for the electronic wave functions, and the 
corresponding distortions of the Fermi surface from a sphere, as well as anisotropy effects 
due to the phonons and umklapp processes in the electron phonon interactions. The Figure 
illustrates the gap Ao(^, 0) at zero temperature, as a function of 6 for three constant arcs. 
Solid angle regions where the Fermi surface of Pb does not exist are indicated by vertical 
solid lines. It is clear that the pure Pb crystal gap is highly anisotropic, varying by about 
20% over the Fermi surface. As described above, impurities will wash out this anisotropy. 
Nevertheless, such anisotropics can be observed in some low temperature properties, like the 
specific heat. For more details the reader is referred to Ref. ||1 1 7|| . 

The other complication we have mentioned is an energy variation in the EDOS, as seems 
to exist in some A15 compounds. If this energy dependence occurs on a scale comparable 
to ujo, then A^(e) cannot be assumed to be constant, and cannot be taken outside of the 
integrals in Eqs. Such EDOS energy dependence is thought to be responsible for 

some of the anomalous properties seen in A15 compounds — their magnetic susceptibility 
and Knight shift [|123|] , and the structural transformation from cubic to tetragonal ||124H126| . 
Several electronic band structure calculations ||127| - p.30|| also find sharp structure in A^(e) at 
the Fermi level. An accurate description of the superconducting state thus requires a proper 
treatment of this structure. This was first undertaken to understand Tc by Horsch and 



Reitschel [118| and independently by Nettel and Thomas [119 . A more general approach to 



understanding the effect of energy dependence in N{e) on Tc was given by Lie and Carbotte 
||120|| , who formulated the functional derivative 6Tc/6N{e); they found that only values of 
A^(e) within 5 to 10 times around the chemical potential have an appreciable effect on the 
value of Tc. More specifically they found that 5Tc/6N{e) is approximately a Lorentzian with 
center at the chemical potential; the function becomes negative only at energies |e — /i| ^SOT^. 

Irradiation damage experiments illustrate some of this dependency. For example, ir- 
radiation of MosGe causes an increase in Tc [|131|| . Washing out gap anisotropy with the 
irradiation cannot possibly account for an increase in Tc, instead, this result finds a natural 



explanation in the fact that the chemical potential for MosGe falls in a valley [|132|| of the 
EDOS, and irradiation smears the EDOS, thus increasing iV(/i), and hence Tc. 

For details on the formulation of Eliashberg theory with an energy dependent A^(e) the 
reader is referred to the work of Pickett | |133| | and Mitrovic and Carbotte ||134|| , and references 
therein. The energy dependent EDOS affects many properties. To illustrate a typical result 
we show in Fig. 6 the effect of an energy dependent EDOS on the current (I)-voltage (V) 
characteristics of a tunneling junction ||134| , |135| . A detailed discussion of tunneling appears 
in Section 3.3.2. The tunneling conductance is proportional to the electron density of states. 



and is denoted by (j{lu) = Rej 



Fig. 6 shows the difference with the BCS 



conductance, a^uj) / aBcsi^) — 1 vs. u — Ao ||134| , |135|| . Fig. 6a (b) is for a peak (valley) in the 



EDOS at the Fermi level. The solid curves include the effect of an energy dependent EDOS, 
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while the dashed curves do not (the EDOS is approximated by a constant value, N{fi)). In 
these examples the electron phonon spectral density obtained for NbaSn ||136|| is used. 

These differences can be highlighted in another way, shown in Fig. 7 ||134|, p.35(1. Here, 



the "effective" electron phonon spectral density, a'^F{Q)eg, is obtained by inverting the solid 
curves in Fig. 6 under the assumption that the EDOS is constant and equal to N{fi). The 
dashed curves give Shen's original a'^F{Q) while the solid curves are the result of (incorrectly) 
inverting the result obtained with an energy dependent EDOS, but not accounting for it in 
the inversion process itself. The actual EDOS used to generate the I-V characteristic is 
shown in the inset for each figure. It contains a peak in Fig. 7a and a valley in Fig. 7b. 
Clearly a peak introduces a negative tail into a'^F{Q)es, which of course is not present in the 
actual a'^F{Q). For other important modifications the reader is referred to the references. 
The rest of this chapter will focus primarily on the 'standard' theory, using Eqs. (|5^|57D 
with Ao(m) = 1 and Ai{m) = 0. 

All of the equations discussed so far have been developed on the imaginary frequency 
axis. Because practitioners in the field at the time were interested in tunneling spectroscopy 
measurements the theory was first developed on the real frequency axis The 
resulting equations are complicated, even for numerical solution. It wasn't until quite a 

137[| , where, for 



number of years later that numerical work returned to the imaginary axis 



thermodynamic properties, the numerical solution was very efficient [p.38| - |141| . The difficulty, 
however, was that imaginary axis solutions are not suitable for dynamical properties. We will 
return to the interplay between imaginary and real frequency axis solutions as we ecounter 
them throughout the chapter. 



3 The Phonons 

3.1 Neutron Scattering 

When dealing with model Hamiltonians, the phonon dispersion relations (before interaction 
with the electrons) are generally given, and simple: they are Einstein modes, or Debye-like 
modes, for example. A noteable exception is the case where the model contains anharmonic 
forces, in which case even the 'non-interacting' phonon spectrum is unknown. 

In the case of real solids, and in particular metals, the situation is much worse. In this 
case the electrons cannot be ignored, though they can be treated in the Born-Oppenheimer 
approximation. Nonetheless the results require parametrization (with input from other ex- 
periments) and are generally not reliable. Pseudopotential methods ||142| , |143| can be applied 



to this problem, again, with limited success. In contrast, the spectacular success of inelastic 
neutron scattering techniques |]89| , pO| to simply measure the phonon dispersion curves in real 



metals effectively eliminates the need to calculate them quantitatively. Various qualitative 
effects, like the impact of electronic screening to the long wavelength ionic plasma mode ||146|| , 



as well as the existence of Kohn anomalies | J147|| , all due to the presence of electrons, are 
understood theoretically. For detailed results, however, Born-von Karman fits to high sym- 
metry phonon dispersions suffice for an excellent description of the low temperature phonon 
properties. At temperatures of order 10 K, the phonons in most conventional superconduc- 
tors are completely determined, and no longer changing with temperature. Hence, as far as 
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understanding (low temperature) superconductivity is concerned, these higher temperature 
measurements are sufficient. 

The measured dispersion curves, cjq (again, branch indices are suppressed), are summa- 
rized in the frequency distribution 

q 

where is the number of ions in the system, and q is a wavevector which ranges over the 
entire First Brillouin Zone (FBZ), (and imphcitly contains the branch index). It should be 
stressed that this procedure is an idealization; in actual fact a set of 'constant q' scans are 



performed (usually along high symmetry directions). A typical result is shown in Fig. 8 
for Pb, for a set of wavevectors along the diagonal in reciprocal space. Note that the neutron 
counts tend to form a peak as a function of energy transfer (to the neutron), hu. In general 
these peaks have a finite width, i.e. broader than the spectrometer resolution; these are due 
to a variety of effects, for example, anharmonic effects. Nonetheless, because the peaks are 
relatively sharp compared to the centroid energy, (i.e. the phonon inverse lifetimes are small 
compared to their energies), these data are usually presented in the form of Fig. 9, as a set 
of dispersion curves. Fig. 9 does obscure, however, the lifetimes of the various phonons, and 
hence the validity of Eq. (0), where infinitely long-lived phonons are assumed throughout 
the Brillouin zone, is called into question. 

Nonetheless, for most of the Brillouin zone the approximation of infinitely long-lived 
excitations is a good one (hence, the name, phonon), and so the spectrum of excitations can 
be constructed according to Eq. ([6^). Such a procedure relies on coherent neutron scattering. 
An alternative is to use incoherent neutron scattering, whereby one measures the spectrum 
more or less directly. This latter procedure has advantages over the former, but also includes 
multiphonon scattering processes, and for non-elemental materials, weighs the contribution 
from each element differently, according to their varying scattering lengths. The result is 
often denoted the 'generalized density of states' (GDOS). A comparison for a Thallium-Lead 



alloy is shown in Fig. 10 [|144| , |145|] . Also shown is the result from tunneling, to be discussed in 



the next subsection. There is clearly good agreement between the various methods. Amongst 
the two neutron scattering techniques, inelastic coherent neutron scattering produces the 
sharpest features, but requires a model (i.e. a Born- von Karman fit) to extract the spectrum 
F(z/) from the dispersion curves measured along high symmetry directions. 



3.2 The Eliashberg Function, a^F{u): Calculations 

First-principle calculations of the electron-phonon spectral function, a'^F{v) require a knowl- 
edge of the electronic wave functions, the phonon spectrum, and the electron-phonon matrix 
elements between two single-electron Bloch states. A fairly comprehensive review is given 
in Ref. [^. For our purposes, we note that, since the phonon spectrum will come from 
experiment, Eq. (pj| ) requires calculation of g-^^^i^'. It is p^,[88[1 

c/k,k'i =< V^k I e^'(k - k') ■ VV I V^k/ > 



h 



(63) 
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where, for this equation we have included the phonon branch index j exphcitly. The Bloch 
state is denoted | V'k >, and e^(k) is the polarization vector for the (jk)th phonon mode. 
The crystal potential is denoted V, and as one might expect, the electron-phonon coupling 
depends on its gradient. 

Tomhnson and Carbotte ||148|| used pseudopotential methods ||149| , |150|| to compute (7k,k'j 

The phonons were taken from experiment [89,^ 
The result is plotted in Fig 



and, from Eq. (16), a^-F(z/), for Pb. 
151| , |152| | through Born - von Karman fits. 



11, along with 

results from tunneling experiments (to be described below). The agreement is qualitatively 
very good; this provides very strong confirmation of the electron phonon mechanism of 
superconductivity. 

Further details of more modern calculations of electron-phonon coupling constants can 
be found in, for example, Refs. and and references therein. Their reliability appears 



to remain an issue, both with the high temperature cuprates, and perhaps less so with 
the fuUeride and more conventional superconductors. The spirit of these calculations is 
somewhat different than the older ones, in that coupling constants are extracted from the 
phonon linewidths, where it is assumed that the phonon broadening is entirely due to the 
electron-ion interaction (and not, say, anharmonic effects). Allen ||153| , |154| derived a formula 
(Fermi's Golden Rule) for the inverse lifetime, 7q(z^), of a phonon with momentum (and 
branch index) q: 



7q = 27ru;q^ |5(k,k' 

k 



/(ek+q - ^J')- /(ek - ^J') 



5(ek+q + hujq - ek). 



(64) 



where again we have suppressed both phonon branch indices and electron band labels. Using 
this equation, in the approximation that the expression [/(ek+q — /i) — /(ek — fi)] /{fi-ujq) is 
replaced by 5(ek — /i) makes it resemble Eq. (plTf ), so that one can write 



E 



1 7. 



3N ^ 2 ^ ^ ^ 



-5{V - tUq) 
I 



where the second line serves to define a q-dependent coupling parameter: 

3 7q 



A. 



(65) 



(66) 



It is through these relations that coupling parameters are often determined. 

It is worth noting at this point that several moments of the function a^F{v) have played 
an important role in characterizing retardation (and strong coupling) effects in supercon- 
ductivity. Foremost amongst these is the mass enhancement parameter. A, already defined 
in Eq. (EDI); in addition, the characteristic phonon frequency, tuin is given by 



cuin = exp 



2 
A7o 



dv In [v) 



(67) 



Further discussion of these calculations can be found in Refs. |TI|,p8 
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3.3 Extraction from Experiment 

Experiments which probe dynamical properties do so as a function of frequency, which is 
a real quantity. However, the Eliashberg equations as formulated in the previous section 
are written on the imaginary frequency axis. To extract information from these equations 
relevant to spectroscopic experiments, one must analytically continue these equations to 
the real frequency axis. Mathematically speaking, this is not a unique procedure; one can 
often imagine several functions whose values on the imaginary axis are equal, and yet differ 
elsewhere in the complex plane (and in particular on the real axis). For example, replacing 
unity by — exp {[3iujm) , in any number of places in the equations does not affect the imaginary 
axis equations, or their solutions, and yet on the real axis the corresponding number of factors 
— exp {(3uj) will appear. 

Physically speaking, however, the Green functions involved have to satisfy certain con- 
ditions; complying with these conditions determines the function uniquely |Q. This allows 
a unique determination of the analytic continuation of the Eliashberg equations on the real 
axis. This procedure will be discussed in the following subsection, followed by subsections on 
experimental spectroscopies, and how they can be used to extract the Eliashberg function. 



3.3.1 The Real- Axis Eliashberg Equations 

We begin with Eqs. - |^ ). To analytically continue Eqs. (^- |^ ) is trivial; one simply 
replaces the imaginary frequency iujm wherever it appears with uj + i5. The i5 remains 
to remind us that we are analytically continuing the function to just above the real axis; 
it is important to specify this since there is a discontinuity in the Green function as one 
crosses the real axis. A simple replacement of iuJm with uj + i5 in Eqs. ( |^ , p5| ) (leaving 
the summations over m') would in general be incorrect. The correct procedure is to first 
perform the Matsubara sum, and then make the replacement. To perform the Matsubara 
sum, however, one has to introduce the spectral representation for the Green functions, G 
and F. These are given by 

^(k, luom) = r du j12Si^ (68) 
F(k,zcu„) = r du^^^, (69) 
where y4(k, cj) is given by Eq. (|2^) and C(k, cj) is given by a similar relation: 

C(k, uj) = --ImF(k, iu + i5). (70) 

TT 



The spectral representation for the phonons is already present in Eqs. (E2|,H3|). Therefore the 



Matsubara sum can be performed straightforwardly (see, for example, Refs. |J3|,^]), and 
the analytical continuation can be done. Upon integrating over momentum (using, as in Eqs. 
( ^^ - |57|) electron-hole symmetry and a constant (and infinite in extent) density of electron 
states), one arrives at the standard real-axis Eliashberg equations 0,|l3|]. These equations are 
much more difficult to solve than the imaginary axis counterparts. They require numerical 
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integration of principal value integrals and square-root singularities, and the various Green 
function components are complex. In contrast the imaginary axis equations are amenable 
to computers (the sums are discrete) and the quantities involved are real. Moreover a 
considerable number of thermodynamic and magnetic properties can be obtained directly 
from the imaginary axis solutions. 

The discrepancy in computational ease between the two formulations led to an alternative 
path to dynamical information, namely the direct analytic continuation of the solutions of 
the imaginary axis equations to the real axis by a fitting procedure with Pade approximants 
| 155| . This method is in general very sensitive to the input data, and has (surmountable 



156| , |157|| ) difficulties at high temperatures and frequencies. 



More recently yet another procedure was formulated ||158|| , which first requires a numer- 
ical solution of the imaginary axis equations, followed by a numerical solution of analytic 
continuation equations. This latter set is formally exact (i.e. no fitting required) and yet 
avoids the complications of the real-axis equations. These equations are 
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f{y - z) + N{y)\ G(k', z - v) ^\f{v ^ z) ^ N{y)\ G(k', z + i/) 1(71) 



kk' 



N 



E 







iV(/i) 



iV(/i) 

'f{u~z) + N{u)]F{k',z 



^ + 



f(u + z) + Niiy) F(k',^ + z/) ,(72) 



where z can actually be anywhere in the upper half-plane. Thus, for example, Eqs. (42,43) 
can be recovered by substituting z = iuJm- On the other hand, once these equations have 
been solved, one can substitute z = u+iS, and iterate the resulting equations to convergence. 
When the "standard" approximations for the momentum dependence are made (i.e. Fermi 
surface averaging, constant density of states, particle-hole symmetry, etc.) the result is 



1 H 2^ A(u; — icUr^ 



UJrr,.Z(iuJr. 



m=—oo 



+ — diya'F{iy){[N{iy) + f{u-iu)] 
(jj Jo 



{uj — u)Z{uj — u + i6) 



+ [N{u) + fiu + co)] 



LU — vYZ'^{uj — u + i5) — (p'^iuj — v + i5) 
{oj + i')Z{ijj + u + iS) 1 



m=—oo 

+m r dva^F{v)\ [N{v) + f{v - u)] 



(73) 



[lUJr. 



(pijjj — u + i5) 



^{uo - uyZ^iuj -iy + i5)- 02(cu -iy + i5) 
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Note that in cases where the square-root is complex, the branch with positive imaginary 
part is to be chosen. 

One important point has been glossed over in these derivations. Because of the infinite 
bandwidth approximation, an unphysical divergence occurs in the term involving the direct 
Coulomb repulsion, Vk,k'5 both in the imaginary axis formulation, Eq. (pBI), and in the 



real-axis formulation, Eq. ([7^). The solution to this difficulty is to introduce a cutoff in 
frequency space (even though the original premise was that the Coulomb repulsion was 
frequency independent) , as is apparent in the two equations. In fact, this cutoff should be of 
order the Fermi energy, or bandwidth. However, this requires a summation (or integration) 
out to huge frequency scales. In fact one can use a scaling argument P, |159| , p!60| to replace 
this summation (or integration) by one which spans a small multiple (~ 6) of the phonon 
frequency range. Hence the magnitude of the Coulomb repulsion is scaled down, and becomes 

m 

where U is a. double Fermi surface average of the direct Coulomb repulsion. This reduction is 
correct physically, in that the retardation due to the phonons should reduce the effectiveness 
of the direct Coulomb repulsion towards breaking up a Cooper pair. It does appear to 



overestimate this reduction, however ||161|| . The analytic continuation of this part of the 
equations has been treated in detail in Ref. ||162|| . 

In the zero temperature limit, Eqs. (^,|7^) are particularly simple. Then the Bose 
function is identically zero and the Fermi function becomes a step function: f{v — uo) 
6{(jj — v). Once the imaginary axis equations have been solved, solution of Eqs. ( [75|J7^ no 
longer requires iteration. One can simply build up the solution by construction from uo = Q 
(assuming a^F{v) has no weight at = 0); in fact, if the phonon spectrum has no weight 
below a frequency, t'min, then only the first lines in Eqs. (^,0) need be evaluated. In 
particular, if the gap (still to be defined) happens to occur below this minimum frequency 
(often a good approximation for a conventional superconductor) then the gap can be obtained 



in this manner 163 



In the following two sections we explore the possibility of using Eqs. ( |75| , [71D to obtain 
information about the microscopic parameters of Eliashberg theory. 



3.3.2 Tunneling 



Perhaps the simplest, most direct probe of the excitations of a solid is through single particle 
tunneling. In this experiment electrons are injected into (or extracted from) a sample, as 
a function of bias voltage, V . The resulting current is proportional to the superconducting 



density of states p8|, |164|-p!66 



Is{V) oc j dojRe 



\uo\ 



CU2- A2(cj) 



[f{uj)-f{u + V)] 



(76) 
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where we have used the gap function, A(ci;), defined as 



A{uj) = (P{u + i6)/Z{Lj + i6). (77) 

The proportionahty constant contains information about the density of states in the electron 
suppher (or acceptor), and the tunnehng matrix element. These are usually assumed to be 
constant. If one takes the zero temperature limit, then the derivative of the current with 
respect to the voltage is simply proportional to the superconducting density of states. 

Re I , 1 , [n 




V2 - A2(V)_ 

where S and N denote "superconducting" and "normal" state, respectively. The right hand 
side of Eq. (|78]) is simply the density of states, computed within the Eliashberg framework 
(see, for example, Ref. [^). It is not at all apparent what the structure of the density of 
states is from Eq. (fTHD, until one has solved for the gap function from Eqs. ( |73| , |7^ and 



Eq. (ffTI). At zero temperature the gap function A{uj) is real and roughly constant up to a 
frequency roughly equal to that constant. This implies that the density of states will have 
a gap, as in BCS theory. At finite temperature the gap function has a small imaginary part 
starting from zero frequency (and, in fact the real part approaches zero at zero frequency 
| 167|| ) so that in principle there is no gap, even for an s-wave order parameter. In practice, a 



very well-defined gap still occurs for moderate coupling, and disappears at finite temperature 
only when the coupling strength is increased significantly ||168| , |169| . 



In Fig. 12 and 13 we show the current- voltage and conductance plots for superconducting 



Pb, taken from McMillan and Rowell [^]. These data were obtained from a superconductor- 
insulator-superconductor (SIS) junction, with Pb being the superconductor on both sides of 
the insulating barrier, so that, rather than directly using Eq. ([78|) , the current is given by a 
convolution of the two superconducting densities of states. Two features immediately stand 
out in these plots. First, a gap is clearly present in Fig. 12, given by 2Ao, where Ao is the 
single electron gap defined by 

Ao = ReA(cu = Ao), (79) 

a definition one can use for all temperatures. Secondly, a significant amount of structure 
occurs beyond the gap region, as is illustrated in Fig. 13. 

McMillan and Rowell were able to deconvolve their measurement, to produce the single 
electron density of states shown in Fig. 14. Since the superconducting density of states is 
given by the right hand side of Eq. (^), the structure in the data must be a refiection of the 
structure present in the gap function, A (a;). The structure in the gap function is in turn a 
refiection of the structure in the input function, a'^F{i'). In other words, Eqs. (^,0) can be 
viewed as as a highly nonlinear transform of a'^F{i>). Thus the structure present in Fig. 14 
contains important information (in coded form) concerning the electron-phonon interaction. 
One has only to "invert" the "transform" to determine a^F(i/) from the tunneling data. 
This is precisely what McMillan and Rowell |^,^ accomplished, first in the case of Pb. 



The procedure to do this is as follows. First a "guess" is made for the entire function, 
a^-F(z/), and the Coulomb pseudopotential parameter, fi*. Then the real axis Eliashberg 
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equations ((|72D and (|73|)) are solved, and the superconducting density of states (Eq. 
is calculated. The result attained will in general differ from the experimentally measured 
function (represented, for example, by Fig. 14); a Newton- Raphson procedure (using func- 
tional derivatives rather than normal derivatives) is used to determine the correction to the 
initial guess for a^F{i/) that will lead to better agreement. Very often another parameter 
(for example, the measured energy gap value) is used to fit /i*. This process is iterated until 
convergence is achieved. The result for Pb is illustrated by the dotted curve in Fig. 11. 

Once a^F{v) (and /i*) has been acquired in this way one can use the Eliashberg equations 
to calculate other properties, for example, T^. These can then be compared to experiment, 
and the agreement in general tends to be fairly good. One may suspect, however, a circular 
argument, since the theory was used to produce the spectrum (from experiment), and now the 
theory is used as a predictive tool, with the same spectrum. There are a number of reasons, 
however, for believing that this procedure has produced meaningful information. First, the 
spectrum attained has come out to be positive definite, as is required physically. Second, the 
spectrum is non-zero precisely in the phonon region, as it should be. Moreover, it agrees very 
well with the calculated spectrum. Thirdly, as already mentioned, various thermodynamic 
properties are calculated with this spectrum, with good agreement with experiment. Finally, 
the density of states itself can be calculated in a frequency regime beyond the phonon region, 
as is shown in Fig. 15. The agreement with experiment is spectacular. 

None of these indicators of success can be taken as definitive proof of the electron-phonon 
interaction. For example, even the excellent agreement with the density of states could be 
understood as a mathematical property of analytic functions ||1 70|| . Also, we have focussed on 



Pb; in other superconductors this procedure has not been so straightforward. For example. 



in Nb a proximity layer is explicitly accounted for in the inversion |p.66| , |171|| , thus introducing 
extra parameters. In the so-called A15 compounds (eg. NbsSn, VaSi, etc.), although the 
measured tunneling results have been inverted ||172|| , several experiments do not fit the overall 
electron-phonon framework [ITD]. 



More details are provided in Ref. [Tl|. An alternate inversion procedure is also provided 
there ||173|| , which utilizes a Kramers- Kronig relation to extract A(ci;) from the tunneling 
result. An inversion of \m.(f){uj + i5) then removes /i* from the procedure. A variant of 
this, where the imaginary axis quantity A{iujm) is extracted directly from the tunneling I-V 
characteristic, and then the imaginary axis equations are inverted for a^F(z/), also works 



174 1, but the accuracy requirements for a unique inversion are very debilitating. 



3.3.3 Optical Conductivity 

In principle, any spectroscopic measurement will contain a signature of a^F(i/). In particular, 
several attempts have been made to infer a'^F[v) from optical conductivity measurements in 



the superconducting state ||175H1771] . In this section we describe a procedure for extracting 
a^F{v) from the normal state ||1 78|| . 

A common method to determine the optical conductivity is to measure the reflectance 
179| as a function of frequency, usually at normal incidence. The reflectance, -R(z/), is defined 



as the absolute ratio squared of reflected over incident electromagnetic wave amplitude. The 
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complex reflectivity is defined by 



r z/. 



= R^/\u)exp{ie{u)), (80) 



where ^(i/) is the phase, and is obtained through a Kramers-Kronig relation from the re- 
flectance ||179| | 

TT JO — V' 

The complex reflectivity is related to the complex index of refraction, n(z/), 

1 — niv) , , 

1 + r;,(z/) 

which, finally, is related to the complex conductivity, oiv) (using the dielectric function, 

e{y) ^ n\v) = e^^ 1^^, (83) 

where is the dielectric function at high frequency (in principle, for infinite frequency this 
would be unity). It is through such transformations that the 'data' is often presented in 'raw' 
form. Nonetheless, assumptions are required to proceed through these steps; for example, 
Eq. (|8lD indicates quite clearly that the reflectance is required over all positive frequencies. 
Thus extrapolation procedures are required at low and high frequencies; a more thorough 
discussion can be found in ||180|| ; see also [ p.81|| . 



For this review, we will consider both static impurities and phonons as sources of electron 
scattering. Both contribute to the optical conductivity, and can be treated theoretically ei- 



ther with the Kubo formalism or with a Boltzmann approach |8^. In the Born approximation 



the result for the conductivity, in the normal state, at zero temperature, is [|176|| : 



V 



where 



cy{v) = ^- du^- r (84) 

roo Q (jj r\i^\ 

T.{uj)= dna^F{n)\ii\- |-ivr/ dVta^F{n) (85) 

Jo il + UJ Jo 

is the effective electron self-energy due to the electron-phonon interaction. The spectral 
function that appears in Eq. (^) is really a closely related function, as has been discussed 
by Allen ||176|| and Scher [p.82|| . For our purposes we will treat them identically. The other 
two parameters that enter these expressions are the electron plasma frequency, tup, and the 
(elastic) electron- impurity scattering rate, 1/r. 

Equation (0) has been written to closely resemble the Drude form, 

the equation could well be recast in this form, with a frequency-dependent scattering rate 
and effective mass (in the plasma frequency) ||183|| . Eqs. (|84D and ( ^5]) make clear that the 
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optical conductivity is given by two integrations over the electron-phonon spectral function. 
One would like to "unravel" this information as much as possible before attempting an 
inversion, so that, in effect, the signal is "enhanced". To this end one can attempt various 
manipulations ||184H186 . 



As a first step one can make a weak coupling type of approximation to obtain ||1 78|] the 
explicit result: 

a^Fiu) = ^^^\uRe^\. (87) 

Note that the conductivity data, including a measurement of the plasma frequency, provides 
us with both the shape and magnitude of a^F(z/). Eq. ( [57| ) works extremely well, as Fig. 16 
shows, in the case of Pb. It tells us that, with a judicious manipulation of the conductivity 
data, the underlying electron-phonon spectral function emerges in closed form. The very 
simple formula, Eq. (^) introduces some errors — it was derived with some approximations 



— as can be seen in Fig. 16. In fact, a full numerical inversion will also succeed ||187| , |188 
the first reference requires a Newton-Raphson iteration technique, while the second uses an 
adaptive method (in the superconducting state). 

Eq. ( P7D was first applied to K^Cqq ||178[| to help determine whether or not this class 
of superconductor was driven by the electron-phonon interaction. The result is shown in 
Fig. 17 and provides convincing evidence that the alkali-doped fullerene superconductors 
are driven by the electron-phonon mechanism. We will return to these superconductors in a 
later section, and further examine the optical conductivity in the superconducting state in 
another section. 



4 The Critical Temperature and the Energy Gap 

Perhaps the most important property of a superconductor is the critical temperature, Tc. 
For this reason a considerable amount of effort has been devoted both towards new materials 
with higher superconducting and, on the theoretical side, towards an analytical solution 
of the linearized Eliashberg equations (set 0^' to zero, where it appears in the denominator 
in Eqs. (|5^ - ) for Tc (see [11, 13| for reviews); the experimental 'holy grail' has enjoyed 
some success, particularly in the last 15 years; the theoretical goal has had limited success. In 
fact numerical solutions are so readily available at present, that the absence of an analytical 
solution is not really debilitating to understanding Tc. 

In the conventional theory there are two input "parameters": a function of frequency, 
a^F(i/), about which we have already said much, and fi*{uJc), a number which summarizes 
the (reduced) Coulomb repulsion experienced by a Cooper electron pair. The focus of this 
chapter will be the effect of size and functional form of a^-F(z/) on Tc. 
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4.1 Approximate Solution: The BCS Limit 

The first insight into Tc comes from reducing the Ehashberg theory to a BCS-hke theory. 
This is accomphshed by approximating the kernel 



2va'^F{v) 



Z/^ + [Ujrn — 



by a constant as long as the magnitude of the two Matsubara frequencies are within a 
frequency rim of the Fermi surface [|140|] , taken for convenience to be Wc, the cutoff used for 
the Coulomb repulsion, /i*. That is, 

X/- • \ /A for both \uJrn\,\^m'\ < ^c, fQa\ 

L U otherwise, 

where A = A(0) has already been defined in Eq. (|20|). Then, the linearized version of Eq. 
(1^) (with Aq^iti') = 1), for the renormalization function, Z{iuJm), reduces to 

Z{iujm) ^ 1 + A. (90) 

Using this and solving the linearized version of Eq. (|56D for the pairing function yields 

41-^4), (91) 



A-/i* ^\27iTc 1) ^^1' 

where '?/'(x) is the digamma function. The cutoff in these equations is along the Matsubara 
frequency axis; this procedure is to be contrasted with the BCS procedure, which introduced 
a cutoff in momentum space. The former is more physical, insofar as the true electron-phonon 
interaction comes from retardation effects, which occur in the temporal domain; hence the 
cutoff should occur in the frequency (either real, or imaginary) domain. In practice, the two 
procedures are connected, so they produce the same physical equation in the weak coupling 
limit. 

Returning to Eq. (|9lD , for large x, il){x) ~ log (x), so, in the weak coupling limit (Tc << 
Wc), we obtain a BCS-like equation, 

Tc = l.lScJcexp f- ^ ^ ^ y (92) 



A - /i* 

This equation has essentially summarized all the detailed information contained in the 
electron-phonon spectral function oP'Fiy) into two parameters, A and uJc- The mass en- 
hancement parameter. A, is a simple moment of o?F{v) (see Eq. (PO])), while the parameter 
LOc physically is meant to represent some typical phonon frequency. In more refined treat- 
ments ||140| , |191[] , LOc is given by some moment of o?F{y) as well. For example, in Ref. ||140|| , 
the logarithmic average is used to define a;in (see Eq. (|67|) ), a quantity we shall use extensively 
in the following sections. They modified the McMillan equation ||191|| to read 

, rj. ^^m ( 1-04(1 + A) ^ 
A derivation of this equation is given in Refs. ||140| , |191| . 
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4.2 Maximum T^, Asymptotic Limits, and Optimal Phonon Spec- 
tra 



Eq. (^) (or Eq. (p3D) describes the weak coupling limit of Eliashberg theory reasonably 
well. It errs in the strong coupling limit; for example, it predicts that Tc saturates as A 
increases, whereas the Eliashberg equations themselves predict that Tc grows indefinitely 
with A [|192|| . Asymptotic results from Eliashberg theory can be obtained correctly and 
analytically [pH], p.4CI| , [193|| through a variety of arguments. The methodology based on scaling 
theorems is particularly powerful, and has been applied to other thermodynamic properties 



as well 



The correct asymptotic result for Tc is ||140| , |193| 



Tc = 0.1S3^/Xue, 



(94) 



as A — s> oo. In obtaining Eq. (p^, an Einstein spectrum has been assumed (this is not 
required), which in turn is characterized by two parameters: the weight, A = Xuje/'^i and the 
frequency, oje- In writing Eq. (plj), one tacitly has assumed that the parameter A = 2A/u}e 
is increased while keeping the frequency uje fixed. In reality, the two parameters are not 
independent — this is the main point of the article by Cohen and Anderson | |194| |. For 
example, often phonon softening occurs because the coupling strength increases. In fact, this 
is made explicit in McMillan's definition |[L91|| of A: 



A 



N{eF)a^ 



(95) 



where N{eF) is the electron density of states at the Fermi energy, M is the ionic mass, and 
is the electron-ion coupling referred to in the Holstein Hamiltonian, Eq. (^) (in a more 
realistic electron-phonon Hamiltonian, would be given by a Fermi surface average of the 



electronic matrix element of the change in crystal potential as one atom is moved 

To determine what the optimal phonon frequencies actually are, functional derivatives 
were introduced 



138 



These had already been utilized extensively as an iterative aid in 
inverting tunneling data with the Eliashberg equations |52|. The most commonly used 
functional derivative is that of Tc with respect to infinitesmal changes in a;^-F(z/), with fixed 



area, A = dv F (u) , first computed by Bergmann and Rainer 
result, derived in Ref. [ p.95|| , is given by the expression 



138|| . An approximate 



6Tr 



E 



4Q 



where Cl = Q/ksTc, and the Bn are numbers given by 

2 2 



(96) 



5. 



Eli 



m=l 



n 2m — 1 (2m — 1)^ 



TT 

T 



(97) 



This function (which is universal) is shown in Fig. 18, and reflects well the generic behaviour 
of the more precise calculation. It illustrates that the optimal phonon frequency lies at some 
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finite frequency (i.e. non-zero, and non-infinite), wliicli is a factor of order 10 times tlie 
critical temperature. Tlius if one could imagine shifting small amounts of weight in a^F(z/) 
around then would increase if spectral weight is shifted either from very high or from very 
low frequencies towards frequencies near the maximum of the curve shown in Fig. 18. 

The reasoning above leads naturally to the concept of an optimum spectrum, first deter- 
mined by Leavens ||196|| , and elaborated upon in Refs. ||197 -201]. In an optimum spectrum 



calculation, one imagines having a fixed area of Q!^F(i/), and asks at what frequency it would 
best be situated in order to optimize some particular property. An appropriate scaling of the 
linearized Eliashberg equations for an Einstein phonon spectrum with frequency ue leads to 
the result 

TJA= f{ujE,l^*), (98) 

where A is the area, ue = uje/A, and / is a universal function of ue, to be determined 
numerically for each choice of fi* (a very weak A dependence in the cutoff associated with 
fi* has been neglected). The result is a curve with a maximum at uje ~ 1; placing a spectral 
function at this frequency will yield the maximum Tc. This procedure yields a result, 

Tc < Acifx*), (99) 

where c{fi*) is a function of n* shown in Fig. 19. Also shown are data from many supercon- 
ductors for which a^F(i/) is known from tunneling spectroscopy, all of which fall below the 
optimum curve. Interestingly, some superconductors have a critical temperature reasonably 
close to their optimal value. 

The last few paragraphs demonstrate the usefulness of functional derivatives in under- 
standing the systematics of Tc. A variant of these results can easily be obtained, which 
may shed even more light on Tc systematics. As we have already mentioned, the functional 
derivative discussed involves the moving around of spectral weight, subject to the condition 
that the area remain constant. However, as Eq. ( pS] ) suggests, it is not the area which likely 
remains constant while phonons soften, but rather the area times a frequency. Hence, one 
can define a different spectral function, 

a^G{u) = ua^F{u), (100) 

and take functional derivatives with respect to this new function. The result is easily obtained 
from that in Eq. (0), simply by dividing by fi. Then ^^2g'(^) will peak at zero frequency, and 
it would seem that it is always advantageous to decrease the phonon frequency. Continuing 



this process will result in a spectrum for which the calculation which gives Eq. (pq ) is 
no longer valid, and one would have to self-consistently calculate the functional derivative, 
numerically. To our knowledge this has not been done for or any other superconducting 
property. 



4.3 Isotope Effect 

As already remarked in the Historical Developments subsection, the discovery of an isotope 
effect on Tc played an important role in the subsequent development of the theory. In the 
BCS equation the isotope effect is clear from the prefactor; phonon frequencies for elemental 
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superconductors are inversely proportional to the square root of the ionic mass, and hence 
the isotope coefficient (3 is 

p^-^ = i (101) 



The last equality follows from Eq. (^), using the fact that A, as defined there, is independent 
of ion mass. In the standard Eliashberg theory, A as defined by Eq. ( PU[ ) remains independent 
of ion mass, and, with fi* = 0, we once again obtain (3=1/2. Complications can arise, for 
example, from a finite electronic bandwidth p02||, or from a non-constant density of states 



near the Fermi level 203-205 



There are two other clear sources of deviation from (3 = 1/2. One is that in non- 
elemental superconductors, an isotopic substitution for one of the elements will result in 
varying changes in T^, depending on how the element being substituted contributes to the 
important phonon modes. One then has to define partial isotope coefficients, defined by 
12061 

^ d InTr , , 

where Mj refers to the mass of the ith element. The total isotope coefficient, Ptot = 

will sum to 1/2 (in the absence of /i*). The total isotope coefficient can also be broken down 

by frequency, with 

and then _ 

Act = / dujf3{uj). (104) 
Jo 

Eq. ( |104| ) is useful when phonon modes coming from one of the elements are well sepa- 
rated from those coming from the others, as exists, for example, in the high temperature 
perovskites, since oxygen is much lighter than the other elements, and hence is chiefly re- 
sponsible for the high frequency modes. Example calculations can be found in Refs. |]lT], |207| . 



The second source of deviation from (3 = 1/2 is because ij*{uJc) is non-zero. To understand 
why this causes less of a reduction in Tc (when a heavier mass is substituted) recall that 
/i*(co'c) is reduced from some larger value nioJB) through the pseudopotential effect. That is, it 
is through retardation that a weaker electron-phonon interaction can overcome the stronger 
direct Coulomb repulsion. In analytical treatments this is often modelled by endowing a 
mass dependency to the Coulomb pseudopotential through the cutoff ||191|| . For example. 



inspection of Eq. ( [75| ) shows a mass dependency if the cutoff frequency ujc is made to 
correspond to a phonon frequency. Then one can derive, from the McMillan equation for Tc, 
Eq. (P3D, the following expression for the isotope coefficient (assuming one element): 

1/ 104(1 + A)(1 + 62A) y 

^ 2V [A-/i*(l + 0.62A)]2 ^ ) ^ ' 

This result properly reduces to 1/2 when /i* = 0, and shows that the isotope coefficient 
is generally reduced when /i* is finite. In fact it is clear from Eq. ( p.05|) that the isotope 
coefficient is reduced for both positive and negative /i*. The isotope coefficient is reduced for 
positive II* because, when you lower the ionic mass, the increase in phonon frequency will 
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raise T^, but not as much as would be the case if the Coulomb repulsion were not present. This 
is because the discrepancy in frequency scales has been reduced slightly, and the retardation- 
induced attractive interaction is not as large as before the isotopic substitution. On the 



other hand a negative /i* represents some unknown attractive non-phonon mechanism ||207 
which contributes to Tc. A weakening of the phonon-induced attraction (through an isotopic 
substitution) reduces only partially, resulting again in a reduced isotope coefficient. 

A more accurate determination of the isotope coefficient can be obtained simply numeri- 
cally, following the prescription of Rainer and Culetto |P06|| . In this case a cutoff is imposed 
on the Eliashberg equations, which is independent of (but much greater than) the maximum 
phonon frequency. An isotopic substitution results in only a shift in the phonon spectrum, 
and a subsequent calculation of will yield the isotope coefficient. This is physically more 
transparent than the analytical approach described above, as an isotope substitution does 
not alter (at this level of theory) the direct Coulomb repulsion. 

There is a substantial literature on the isotope effect; much of the older results are 
summarized in Ref. |[164|| . By 1969 many low Tc superconductors had been found, several 



of which had very low isotope coefficients. These had, for the most part, been explained 
through detailed calculations P, |208|| , due to the physics outlined above. It is worth noting 



that this explanation of the sometimes low isotope coefficient observed was not universally 
accepted ||209| . 



The discovery of the high temperature cuprate materials prompted considerable activity 



concerning the isotope coefficient, as is reviewed in Ref. [piCI|| . The isotope coefficient displays 
some unusal doping dependence in the La2-x{Sr,Ba}a;Cu04_5 systems, but is essentially zero 
in the optimally doped 90 K YBa2Cu307_y system. The question is, can a realistic (and 
conventional) electron-phonon interaction give rise to a 90 K superconductor with a near 
zero isotope coefficient ? A qualitative answer can be obtained ||211|| through the use of the 
McMillan equation ( p^ , |105|) . For a given electron-phonon coupling. A, and phonon frequency 
uje, one can determine the required value of fi*{ujc) to fix Tc from Eq. ( PB| ) (assuming Uc 
refers to a cutoff associated with the phonon spectrum). These parameters can then be used 
in Eq. ( |105|) to determine the isotope coefficient, /?. Some such results are plotted in Fig. 



20. To obtain the desired results for optimally doped YBa2Cu307_y (Tc 90 K and /3 ~ 0) 
would require high frequency phonons uje ~ 100 meV with very strong electron-phonon 
coupling (A ~ 5). That such a coupling strength is unrealistic, particularly for such very 



high frequency phonons, was discussed much earlier by Cohen and Anderson ||194 . 

The qualitative validity of Fig. 20 has been verified by several numerical solutions to the 
Eliashberg equations ||202| , |21l| , pl2|| . In particular, in Ref. ||202|| a natural bandwidth cutoff 



was employed, with similar results. In summary the conventional Eliashberg theory can yield 
a near- zero isotope coefficient, provided Tc is low. One must go beyond the conventional 
framework to obtain a zero isotope coefficient with Tc ~ 90 K. 

4.4 The Energy Gap 

The existence of a single particle energy gap, although not fundamental to superconductivity 
|213|, nonetheless has played an important role in our understanding of superconductivity. 



How an energy gap arises in the I-V characteristic of a conventional superconductor has 
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already been discussed in Section 3.3.2; there we focussed on extracting detailed information 
about the mechanism. Here we turn our attention to the gap, a much more prominent 
feature in the experimental result, and learn what a particular value may imply about the 
superconductor. 

The first step is to examine what occurs in BCS theory. The order parameter is then 
given by a constant, as written in Eq. (^). Suitable generalization to the model interaction 
given by Eq. ( |gPD yields 

A = 2cu,exp(^-^^), (106) 

in the weak coupling limit. The solution at finite temperature is somewhat more complicated; 
it can be obtained numerically, and shows the typical mean field behaviour near Tc [§]: 

A(T) oc (Te - r)^/^ (107) 

Near T = the order parameter is exponentially fiat ||214|| : 



A(T) ^ A(0) - [27rA(0)T]i/2g-A{o)/r_ ^^Qg^ 



The order parameter is a real (i.e. not complex) number for all temperatures ||215|| . Thus 



Eq. ([78|) shows that the dl/dV curve (which provides an image of the density of states) 
will show an energy gap at A(T) at each temperature. An illustration of the temperature 
dependence of the order parameter is given in Fig. (21a), along with the density of states at 
several temperatures (Fig. (21b)). 

Within Eliashberg theory, the calculation of the corresponding property is much more 
complicated. First of all, a careful distinction between the gap or pairing function (which 
is now a function of frequency at any given temperature) and the energy gap is required. 
The energy gap is defined through Eq. (^). The gap function is, in general, a comphcated 
and complex function of frequency, that results from a solution of the Eliashberg equations. 
These, in turn, can be solved either on the imaginary axis (Eqs. ( ^If) and (^) for the 
'standard' theory) or the real axis (Eqs. (^) and ([74D). Example solutions for a real 
electron-phonon spectrum (Pb) are shown in Fig. 22 and Fig. 23. The solutions on the 
imaginary axis turn out to be real; on the real axis they are complex. The corresponding 
densities of states are shown in Fig. 24. 

The low frequency behaviour of the various functions plotted is not clear on the figures 
shown. A careful analysis ||167|| leads to 



Re A (a;) = c 
lmA{uj) = 
ReZ^u) = d 
ImZioj) = 



T = 



T = (109) 



at zero temperature, where c and d are constants, whereas at any non-zero temperature, we 
obtain 

ReA(w) oc ^2 ^ ^ ^ 
ImA(co') oc u 

ReZ(uj) = d(T) , , 

T \ , T>0. 110 

ImZ oc 1/cj ^ ' 
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The latter result in particular implies that, strictly speaking, at finite temperature there 
is always "gapless" superconductivity. However, as can be seen from Fig. 24, in reality 
the "finite temperature" density of states at zero energy is generally quite small (except for 
very close to Tc). The extent to which this is true depends on the electron-phonon coupling 
strength; as this increases the zero frequency density of states can be a significant fraction 
of the normal state value at temperatures near [|169|| . 

In the remaining subsections we wish to examine the dependence of the energy gap 
on coupling strength. Since the electron-phonon interaction is characterized by a spectral 
function, a^F(fi), we first must decide how to quantify the coupling strength of a particular 
superconductor. Historically the mass enhancement parameter, given by Eq. (pOD, has played 
this role. However, depending on the material, the direct Coulomb repulsion, characterized 
by fi*{ujc), where Uc is some suitable cutoff frequency, can offset the effect of A. Another 
possible parameter is the ratio of the critical temperature to an average phonon frequency. 



a quantity first advocated by Geilikman and Kresin [21f:-219|. This approach was further 



quantified by Mitrovic et al. ||220|| . In this reference (see also Ref. ||221|| ), the Allen-Dynes 
parameter Tc/uj\^ emerged naturally in the derivation of strong coupling corrections, as an 
indicator of coupling strength. A large number of superconducting properties were obtained 
in this way (see Refs. |p!T|, |222|| for derivations and more details), and semi-empirical fits 



were obtained based on accurate numerical solutions. We discuss these further in the next 
section. 



4.5 The Energy Gap: Dependence on Coupling Strength Tc/uj^n 

As we have already emphasized, cannot be reliably calculated at present. The first, 
perhaps simplest, test for the accuracy of Eliashberg theory is then its ability to properly 
obtain the gap ratio, 2Ao//cbTc, where, by Ao, we mean the zero temperature gap edge. In 
Ref. [E2D1 (see also Ref. |TT|), numerically calculated results were compared to experimental 



tunneling results for Ao, obtained for a variety of conventional elemental and alloy supercon- 
ductors. The deviations of the gap ratio from the BCS universal result, H^oI^b^c = 3.53, 
are up to 50 %; yet the level of error is about 5 % , with one notable exception (NbaSn). The 
theoretical results are obtained from a solution of the imaginary axis equations (Eqs. ( ^4|) 
and (|56[), with the standard approximations of infinite bandwidth and particle-hole symme- 
try), followed by an analytical continuation to the real axis. To obtain the gap edge, a Fade 
approximant suffices to get very accurate results ||155|| , as the more systematic continuation 



|158|| verifies. In any event it is desirable to have an analytic form for these corrections. The 



result of Mitrovic et al. [ ^2CI[ | is 

2A„ 



3.53 



l + 12.5(^)^ln(^: 



In obtaining this result the spirit of the McMillan equation was followed, and the coefficients 
12.5 and 2 were chosen from fits to the numerical data for a large number of superconductors. 
These results are plotted in Fig. 25. From this Figure it is clear that Eq. ( p.ll|) describes 
the overall trend very well. As the electron-phonon interaction increases (i.e. becomes more 
retarded), the gap ratio increases to values exceeding 5.0. 
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Figure 25 illustrates that a simple analytic form describes the trend of the gap ratio 
as a function of T(./uj\n rather well for a variety of conventional superconductors. In each 
case electron-phonon spectral functions were used, as obtained from tunneling data, or, 
in some cases, model calculations. On occasion, one sometimes uses a phonon spectrum 
obtained from inelastic neutron scattering, scaled to give the measured critical temperature. 
This latter process assumes that the electron-phonon coupling is constant as a function of 
frequency (seen to be reasonable in the case of Pb), and often assumes a value of the Coulomb 
pseudopotential, ^*{u)c = ScUmax) ~ 0.1 (cUmax is the maximum phonon frequency). Specific 
references to the sources of these spectra can be found in Refs. |TT| , |222| |. 

An important question, particularly when faced with a new superconductor whose phonon 
characteristics may or may not be 'typical', is to what extent the trend modelled by the semi- 
empirical analytic form, Eq. ( |111| ) can be violated, for a given coupling study. This question 
was considered in Ref. |p23|| . They took existing electron-phonon spectra, alF{v), and scaled 
them to new spectra, alF{v)* = BalF{hv), where B and h are constants, chosen to span 
a continuum of values of Tc/tuin. Thus, given some spectral shape, say that of Pb, one can 
determine a curve of I/S^oI^bFc vs. Tc/cuin- In this way they were able to ascertain, for a given 
value of Tc/tuin, the shape dependence of the gap ratio. They of course found more significant 



deviations from the analytical form, Eq. (|111D; nonetheless, the deviations remained small 



on the scale of Fig. 25. Larger deviations were obtained with the use of (somewhat artificial) 
delta-function model spectra [p23H225| . Similarly, if the electron-phonon coupling strength 



is taken to be extremely high, large deviations occur from one spectral shape to another 
226| . 



The net conclusion is that, with physical spectra and physically relevant coupling strength 
(Tc/co'in^O.2), the strong coupling corrections are quasi- universal, and are well described by 
Eq. ( |111| ). We explore in the next subsection how this can be used to optimize the gap and 



gap ratio. 

4.6 Optimal Phonon Spectra and Asymptotic Limits 

A functional derivative analysis similar to that described for Tc yields, for Ao, an optimum 
phonon frequency for a given spectral area. One finds that for a delta function spectral 
function, the zero temperature gap edge obeys a scaling relation just like given by Eq. 

® AoM = (7(cus,/i*), (112) 



where all quantities are as defined following Eq. (p8|). As found there, for a given base spec- 
trum, an optimum frequency U:*^ exists whose value is generally lower than the characteristic 
frequency of the base spectrum — this is particularly clear when the base spectrum itself 
is a delta function. With one found that shifting the spectral weight to that optimum 
frequency resulted in an enhancement of T^. Furthermore, an iteration of this procedure 
resulted in convergence to the situation where, for a given spectral area, the maximum 
had been achieved, with a frequency given by a;^ ^ 1.3 (for \i*{uJ() = 0.1). The functional 
derivative of Tc with respect to q;^F(z/) using this base spectrum is non-positive definite ||199 
with a maximum at cj^, showing that Tc could no longer be increased. 



c 
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The situation with the gap edge is similar, but differs in the following crucial point. 
Upon iteration one finds that the optimum frequency continues to decrease, as the Einstein 
frequency of the base spectrum decreases. Thus, the implication is that the gap edge, 
and therefore the gapedge ratio, 2Ao/fcBTc, will be maximized in the limit as a;^ 0. 
Alternatively, since these calculations are for fixed spectral area. A, this will occur as A ^ oo. 

What is the maximum value of 2Ao/ ksT^ allowed within 'standard' Eliashberg theory ? 
Carbotte et al. | p,98| | answered this question through a scaling theorem, and backed up with 
numerical work. They found that the gap ratio increased monotonically as A increased, 
finding (numerically) a value close to 10 (recall BCS gives 3.53) for values of A ~ 30. In 
doing so they proved that Aq oc \/\uje as A oo, just like does (Eq. (0)). Claims were 
made to the contrary, but these were definitively put to rest in Ref. ||169|| . By solving a set 
of Eliashberg equations written specifically for A — oo, they found a maximum value of the 
gap ratio equal to 12.7. A variety of other properties were explored in the asymptotic limit, 
A — *^ oo, as can be found in the previous references and in Ref. |p271 - |231 |. 



5 Thermodynamics and Critical Magnetic Fields 



These topics have been amply covered in previous reviews [irT| . Nonetheless, we include here 



for completeness a brief summary of the impact of the electron-phonon interaction on these 
properties in the superconducting state. 



5.1 The Specific Heat 

To calculate the specific heat one requires the free energy. For an interacting electron system, 
a practical formulation of this problem was first proposed by Luttinger and Ward ||232|| , and 
further pursued by Eliashberg 



|. A simpler calculation requires the free energy difference 
between the superconducting and normal state, for which an expression due to Bardeen and 
Stephen ||234|| is 



AF 

iV(0) 



A2| 



- z 



lUJ.n 



(113) 



where, for clarity, we include the two Eliashberg equations from Eqs. (|5^ - 



1 + ttT^ A(iu;m - 



TxT^[\{iuJm - iuJni') - A^(0)Vcoul 



(114) 
(115) 



These equations ignore band structure effects entirely (except through the electron density of 
states at the Fermi level, denoted here by A^(0)), and we again have adopted the shorthand 
Z{iu}m) = Zm etc., and used the gap function A{iuJm) = 4>{ii-0m) / Z {iujm) ■ For the free energy 
expression we have used superscripts 'S' or 'N' to denote the superconducting or normal 



39 



state, respectively. In the normal state Z^iiujm) reduces to the expressions obtained in 
subsection (2.3), which is easily seen if one uses the relation 

nz)=z{l-Z{z)), (116) 

where 2 is a frequency anywhere in the upper half plane, and ^(2) is the electron self-energy. 
Equation ( |113| ) can easily be evaluated, once the imaginary axis Eliashberg equations 



(114- 115 ) are solved. From this the specific heat difference. 



j2 A p 

AC(T) = -T^^, (117) 



and the thermodynamical critical field. 



H,{T) = V-SttAF, (118) 

can be computed. The former displays a jump at Tc, characteristic of a mean field theory, 
which is the level of approximation of Eliashberg theory. At low temperatures the specific 
heat in the superconducting state should be exponentially suppressed. This is generally ob- 
served [[141|| , and deviations that do occur at very low temperatures can readily be explained 



by anisotropy in the gap parameter [|117]| . 

Because properties like the electron density of states at the Fermi level are difficult to 
measure or calculate reliably, one would like to focus on observables that are independent 
of these properties. For the specific heat difference, one way of accomplishing this is to 
normalize the specific heat to the normal state result, which presumably contains the same 
electron density of states. The result is then independent of A^(0), and can be compared 
directly to the measured results. A textbook example was provided in the case of Al ||235|| ; 
the data is reproduced in Fig. 26, along with the BCS prediction. The normal state specific 
heat for a weakly interacting electron gas is given by 

C^(T) = 7T, (119) 

where 7 is the Sommerfeld constant given by 

j = ^n'klN{0){l + \). (120) 

Here, A is the electron-phonon enhancement parameter, already referred to on many occa- 
sions. The electron-phonon interaction alters the low temperature specific heat through the 
mass enhancement parameter, 1 + A. In fact, a more careful treatment jSSl , |100| , |236|| yields 



a temperature-dependent 7(T) for the specific heat coefficient (which, at very low tem- 
perature, reduces to the Sommerfeld 7). Besides providing quantitative corrections to the 
electronic specific heat in the normal state, this correction also provides a properly physical 
contribution from the low frequency phonon modes, as found in Ref. [p37|| . 

For a variety of conventional superconductors, like Al, the normal state low temperature 
specific heat is easily measured by suppressing the superconducting state with a magnetic 
field. Then the ratio AC{T)/'~fTc can be determined. At the BCS result for this ratio is 
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universal, like the gap ratio: it is 1.43. Strong coupling corrections can be derived | 221 ] as 
before, as a function of the strong coupling parameter, Tcjuji^. The result is 



7Tc 



1.43 



l + 53(^fln(|;^) 

LO\n Ale 



(121) 



Again, the coefficients 53 and 3 were determined semi-empirically by fits to numerical data. 
A plot of this result, along with some of the numerical data, is shown in Fig. 27. We already 
remarked about Al — its calculated value is indicated by the point nearest the ordinate, and 
agrees very well with experiment. The result for Pb is also shown; the experimental value 
is 2.65, almost a factor of 2 greater than the BCS result. The theoretical result, based on a 
numerical solution of Eqs. (|113[ - [L15| ), is in good agreement. 

In particular, the asymp- 



The result for stronger coupling has also been calculated |226 



totic limit can be computed, following standard procedures. The result is |229 



19.9 



(122) 



showing that the relative magnitude of the jump decreases for large A, and therefore, as is 
already becoming apparent in Fig. 27, the specific heat will have a maximum as a function 
of coupling strength. 

Similar results can be derived for other thermodynamic properties as well. These have 
been summarized in Ref. 



Ill and will be omitted here. 



5.2 Critical Magnetic Fields 

In a type-I superconductor, a critical magnetic field {He) exists, given by Eq. ( |118| ). In a 
type-II superconductor, a lower critical field, ifci, and an upper critical field, ifc2, exist; the 
former signals the departure from the Meissner state to one in which one vortex penetrates 
the system, while the latter occurs at the normal/superconducting transition. The thermo- 
dynamic critical field continues to exist as a thermodynamic property, but not one that can 
be measured by application of a magnetic field. 

A theory of Hd has been worked out within the BCS approximation in Ref. ||238| , p39 | 



(in the dirty limit). This work was extended to the level of Eliashberg theory in Ref. ||240|| . 
It is traditional to calculate the reduced field, hdiT /T^ = y^^*-^^^ as a function of T/Tc. 
Such a curve has a slope of —1 near T^, and saturates to some value at T = 0. Rammer ||240|| 
found that the low temperature value decreased with coupling strength (there characterized 
by a particular spectrum). 

A detailed theory has also been provided for the upper critical field, Hc2- In 1957 
Abrikosov essentially created the subject of type II superconductivity ||241|| . Both exper- 



imental and theoretical work in this exciting area continued to fiourish throughout the 
1960 's. Applications of superconductivity in the mixed state require type II superconductiv- 
ity in order to sustain high magnetic fields. Abrikosov's solution used the phenomenology 
of Ginzburg-Landau theory |^|. Further theoretical developments utilized the microscopic 



theory of Gor'kov [|^. The first of these was by Gor'kov ||242|| for clean superconductors. 
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followed by five papers by Werthamer and collaborators to include impurity effects | 243 , 244 
spin and spin-orbit effects [p45|| , Fermi surface anisotropy effects ||246|| , and retardation effects 
|247|| . All of these papers used an instantaneous attractive potential (i.e. as in BCS theory), 
except for the last. Further developments to include retardation effects were carried out 
in Refs. Q |238|, |48[ and others. Finally, in Ref. [g the Eliashberg theory of i/c2, 
including Pauli paramganetic limiting and arbitrary impurity scattering, was formulated and 
solved. 

Without retardation effects or Pauli limiting, the zero temperature upper critical field, 
when expressed in terms of the slope near T^, takes on universal values, dependent only 
on the elastic impurity scattering rate, given by 1/r. For example, the quantity hc2{0) = 
H^2{^)/{%\H'^^{T^)\) is given by 0.693 in the dirty limit (1/r » A) and 0.727 in the clean 
limit (1/t = 0). For intermediate scattering rates the result falls somewhere in between. It is 
worth mentioning that the absolute value of the upper critical field increases with increased 
impurity scattering. We often use the ratio because the slope near Tc is measured, and 
then the zero temperature value is obtained by using the universal number quoted above. 
The value at zero temperature is of special interest because the Ginzburg-Landau coherence 
length can then be extracted through 



H, 



c2 



(123) 



Here $o is the fiuxoid quantum, and we have used the subscript 'GL' to denote the Ginzburg- 
Landau coherence length, which, at zero temperature, is often close to the BCS coherence 
length, and gives us an indication of the Cooper pair size |P50||. Hence, deviations from 0.693 



(or 0.727) due to retardation effects are of interest for this reason. 

For completeness, we quote the equations which govern Hc2-, taking into account electron- 
phonon interactions in the Eliashberg sense, and Pauli limiting. The gap equation is linear 
in the order parameter [ p49| 



X-H^(^^m))-l/2r' 



with 



m 

The factor xi'^{i^n)) is given by 
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(124) 



(125) 



(126) 



Here a{T) = ^\e\Hc2{T)vp, with e the charge of the electron and vp the electron Fermi 
velocity, /x^ is the Bohr magneton. Eq. ( |124| ) can be written as an eigenvalue equation, just 
like Tc. It is linear because the solution is valid only on the phase boundary between the 
normal and the superconducting states. 

We have carried out extensive numerical investigations of /ic2(0) as a function of coupling 
strength. In the conventional regime, the dependence on coupling strength is very weak 
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251|| ; in the dirty limit hc2{0) decreases initially (as a function of Tc/cuin) to about 0.65, and 



then increases to beyond 0.70. In the clean limit there is first a barely discernible decrease, 
followed by an increase to values of approximately 0.80. These are all theoretical results, 
and, in many cases have not been carefully investigated with experiment. On the other hand 
the expected changes are of order 10% or less, and may well be masked by other effects. A 



thorough investigation was provided for Nb by Schachinger et al. ||252|] . The agreement with 
the available data was excellent (although they did invoke, in addition to the theoretical 
framework described here, anisotropy effects). For further information the interested reader 
is directed to the aforementioned references. 

Before leaving this section we should also mention that optimum spectrum analysis ||201 



and asymptotic limits (3rd reference in [229]) have also been investigated for Hc2] the result 
is very dependent on elastic impurity content, except in the asymptotic limit. In that case, 
the results approach a universal value, i.e. ^c2(0) — >■ 0.57 as A increases. 



6 Response Functions 

In the previous sections we have seen effects due to the inclusion (through the Eliashberg 
formalism) of the detailed electron-phonon coupling. The result is in many cases a large 
quantitative correction to the corresponding BCS result. In this way one can infer, from 
experiment, the necessity of taking into account the dynamics of the electron-phonon in- 
teraction. Nonetheless, as we saw in Section 3 (particularly in the Tunneling and Optical 
Conductivity subsections) dynamical interactions manifest themselves more clearly in dy- 
namical properties. For this reason we now focus on various response functions. 



6.1 Formalities 



A theory of linear response can be approached from two very different frameworks, the 
Kubo formalism, and the Boltzmann equation. The two frameworks often lead to the same 
result; their connection is discussed at length in Ref. |8^. Early treatments P, |253|| of the 
various response functions in a superconductor neglected the electron-phonon interaction, 
except insofar as it provided the mechanism for the superconductivity in the first place. The 
main interest was the investigation of a new state which apparently had a single electron 
energy gap, which would manifest itself either directly in spectroscopic methods (optical and 
tunneling) or more indirectly as a function of temperature (NMR relaxation rate, acoustic 
attenuation, etc.). Sometime later two seminal papers appeared [ |10U| , |254j |, both of which 
discussed the impact of the electron-phonon interaction on transport in the electron gas. 
These dealt specifically with the normal state. Work at a similar level but in the super- 
conducting state appeared a little later [p55|| ; this latter work was generahzed to apply for 
arbitrary elastic impurity scattering only much later ||256|| . These authors used quasiclassical 
techniques; below we will sketch an alternative derivation based on the Kubo ||257|| formula. 

We should preface this work with some remarks about vertex corrections. They are gener- 
ally ignored in calculations of response functions, so that a particle- hole 'bubble', consisting 
of one single electron Green function and one single hole Green function, requires evaluation 
258|| . Older work | ]263| , |264j] investigated the need for vertex corrections and found that they 
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contributed very little; later work in the normal state ||1 76[ p.82|| suggested that their contri 



bution could be summarized by substituting a 'transport' electron-phonon spectral function, 
q;^j.F(z/), for the usual spectral function, a^F(z/), in the transport equations. This alteration 
discriminated in favour of back scattering as being particularly effective in depleting the cur- 
rent, as one would expect. Over a large frequency range, however, these spectral functions 
are not expected to differ substantially; nonetheless, quantitative investigations are currently 
lacking particularly in the superconducting state. 

The contribution to the conductivity consists of two components: the paramagnetic and 
diamagnetic responses. The diamagnetic response is straightforward ||265|| ; the paramagnetic 



response is determined by the evaluation of a current- current response function. A standard 
decoupling of this function (ignoring, as noted above, vertex corrections) yields 



2 

(U{u + i5) + —), (127) 



where n(z/-|-i5) is the paramagnetic response function whose frequency dependence (on the 
imaginary frequency axis) is given by 

n(iz/„) = ^r{ev^yG{k, iUm)G{k, iui^ + iz/„), (128) 

where G{]<.,iu}m) is actually a matrix in the Nambu formalism 0. It is given, in terms of 
functions with which we are already familiar, by 

Lr[ii,lUJjn) — -J T-. 7^. xxo , lO^' V' y^"^^) 

(Ck - yU)^ - {tUJ,nZ{lUJm)r + (P [^^m) 

where the Pauli spin matrices are given by 

?)' ^^"(? J)' ^^"(? 7)' ^^"(J -l)- ^^2°) 



In Eq. ( [1281) the trace is over the Pauli spin space. The presence of the factor {eVxY shows 



explicitly that the dressed vertex has been replaced with a bare vertex; is the component of 
the electron velocity in the x-direction. The momentum sum is over the entire Brillouin zone; 
the factors preceding the summations include the total number of atoms in the crystal, A^, 
and the inverse temperature, (3 = l/ksT. The diamagnetic piece in Eq. (|127|) contains the 
electron density n and the electron mass, m. Finally, the single electron energy is denoted by 
ek, and, as before, we use a notation where we explicitly subtract off the chemical potential. 

The paramagnetic kernel denoted by n(zz/„) in Eq. ( p.28| ) is a special case of a more general 
'bubble' diagram. A similar calculation, for example, is required for the NMR relaxation 
rate |[L68| , |266| , |267|| , or the phonon self-energy |P68H271|| , except that the vertices are not 



proportional to tq (as was the case in Eq. ( |128|) ), but to some other Pauli matrix. This 
has the effect that the so-called 'coherence factors' will differ, depending on the particular 
response function; some will result in a cancellation with singularities arising from the single 
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electron density of states, whereas others will result in a potentially singular response, at low 
frequencies, as we shall see below. An early review outlining these differences in the context 
of Eliashberg theory is given in Ref. [0| . 

Returning to Eq. (|128|) , the standard procedure is as follows; one would like to evaluate 
the Matsubara sum — only then can one perform the proper analytic continuation to real 
frequencies required for the optical conductivity. This is straightforward, through the spec- 
tral representation, which is the Nambu generalization of Eq. The cost is that two 
new frequency integrals are required, one of which can be done immediately by making use 
of the Kramers-Kronig-like relation (see Eq. 



G(k,z)= r dco^^^, (131) 

J —CO Z — UO 

with z anywhere in the upper half plane. Finally, we would like to perform the Brillouin 
zone integration analytically; to do so, we note that the only dependence on wavevector k 
in Eq. (|128|) occurs through ek (this is not so for more complicated response functions, such 
as is required for neutron scattering, for example, where the momentum dependent kernel, 
n(q, u + i6), is required). This feature of the optical response allows us to make the usual 
replacement, given already by Eq. 



l^E^ I deNie), (132) 

k 

where A^(e) is the single electron density of states. As in that case A^(e) can be taken as 
constant (= N{ep)) and, along with the electron velocity, v^, taken out of the integration as 

2 

an overall constant, 2N{eF)e'^vl = ^ = ne'^/m, where up is the electron plasma frequency. 
However, one would normally like to extend the integration over single electron energy from 
— oo to +00, as is often done within Eliashberg theory. Here, however, one has to be slightly 
more careful, and first subtract the normal state contribution to the kernel. This makes 
the integral sufficiently convergent that extension to an infinite bandwidth (effectively) is 
possible. Then the integral can be readily performed by contour integration. The integration 
over the normal state contribution alone must be done separately; an integration cutoff ±D 
must be used, the effect of which is an additional (imaginary) contribution. The final result 

vm 
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N{uj) 
Piuj) 



+ iS) 
e{u! + i6) 



e{uj + id) 

e{uj) = ^u'^{uj + i5)-(l)^{uj + i5) (134) 

where D is the large cutoff mentioned above, to be taken to infinity for large electronic 
bandwidth, and Cj{ijJ + i5) = ujZ{ijj + i5). 

Various limits can be extracted from these expressions; for example the normal state 
results of Section 3.3.3 can be readily obtained, as well as the simple Drude result, obtained by 
assuming only elastic scattering characterized by a frequency independent rate, 1/r. When 
inelastic scattering is included (here through electron-phonon scattering), low frequency 



Drude-like fits can be obtained through simple expansions ||184|| . We will turn to these later. 

Equation ( |133| ) represents the 'standard' theory of the optical conductivity with Eliash- 
berg theory. As already mentioned, this characterization includes the caveats discussed above 
about vertex corrections and a^F{v) —>■ a^j.F(z/) replacements. It is valid for both inelastic 
scattering and elastic scattering processes (within the Born approximation). The impact of 
elastic scattering on the Eliashberg equations have not yet been discussed, so we turn to 
these now. Equations (|114|J115|) , on the imaginary axis, along with Eqs. (^,^), on the 



real frequency axis, are written for the clean limit ('clean limit' is here defined to mean that 
the elastic scattering rate is zero, 1/r = 0). When elastic scattering is included, new terms 
appear on the right hand side of these equations. (As an aside, one way of using the existing 
equations to include elastic scattering is to include a component of q;^F(z/) (called af^pF^u) 
for simplicity) which models the elastic scattering part. At any non-zero temperature it will 
be given by a zero frequency contribution 



Substitution of this expression into Eqs. (|114| , |115D , for example, will yield simple expres 



sions on the right hand side proportional to 1/r.) In principle one would think that Eqs. 
( |114|J115| , [73|j7^) require iteration to a solution for every new value of impurity scattering. 



In actual fact, however, they need be solved only in the clean limit. Then, the pairing 
+ i6), and renormalization u)(a; + i6) = ujZ{uj + i6) functions can be modified by the 
simple contribution 

uico + tS) ^ Coico + 15) + — , ^(^ + ^'^) (137) 

2r ^Cj\uj + i5) -(t)^{uj + i5) 



Equations ( |133| , p.34] ) remain the same with impurity scattering. This is a consequence of 



the so-called Anderson's 'theorem' [|121|| . The modifications are all implicitly contained in 



the pairing and renormalization functions. Note that the gap parameter, A(cij + i5) = 
(j){uj + iS) / Z{uj + ib\ remains the same, independent of the impurity scattering rate. 
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6.2 BCS results 



The purpose of this chapter is to examine effects specifically due to the electron-phonon 
interaction. Nonetheless, it is best to first see what occurs in the BCS limit, and then 
examine the differences. The means for achieving the BCS limit from Eliashberg theory was 
examined in Section 4.1; in general we mean by the 'BCS limit' that limit which corresponds 
to taking a^-F(i/) to be non-zero only for some very high frequency component (so that 
the "strong coupling" indicator Tc/win ^0). As a result, the renormalization function, 
Z{u + i6) 1, and the gap function A{u + i6) ^ A, a constant, as a function of frequency. 
This allows one to explicitly break up the integrals in Eq. ( |133|) into portions involving the 



BCS gap parameter. A, and the electromagnetic frequency, u. A very efficient FORTRAN 



program has been provided in Ref. ||273|| in this case 



6.2.1 Far-Infrared: Dirty Limit 

A historically important case is the dirty limit. This is defined by 1/r >> A, and was 
first treated by Mattis and Bardeen [ ^53| | . An analytical expression can be obtained at zero 
temperature [|253|| : 



— = {l + —)E{k)-—K{k) z/>2A (138) 
^ = hl+^A^E(k-)-l(l-?^)K(k'), (139) 

2 

where (j„ = ^^^^^ is the normal state conductivity (pure real) and the real part of the con- 
ductivity is identically zero for frequencies, v < 2A. In these expressions 

= I ~ I and k' = Vl- (140) 

and E{k) and K{k) are the complete elliptic integrals of the first and second kind. For 
other cases (finite temperature and/or lower impurity scattering rate) one must integrate 
numerically ||273|| . Figure 28 shows (a) the real part and (b) the imaginary part of the con- 



ductivity in the zero temperature BCS superconducting state, for various impurity scattering 
rates. We have used some definite values for the impurity scattering rates and the coupling 
strength. The latter has been chosen to yield an absorption edge, 2A = 10.4 meV, which, 
because of the insensitivity of superconductivity to elastic impurity scattering [ |121|| , holds 
for all scattering rates. A well-defined absorption onset is evident in Fig. (28a); otherwise 
the curves simply deviate from what would have been Drude-like curves in the normal state. 
In Fig. (28b) the frequency times the imaginary part of the conductivity is shown for the 
same scattering rates. Such a combination is shown because the zero frequency limit gives a 
direct measure of the London penetration depth: 

4-Tr 

l/A2(T) = lim— z/a2(z/). (141) 

u^O c 

As is evident from the figure, the penetration depth increases as the impurity scattering rate 
increases. 
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Another feature stands out in Fig. (28b); there is a notable 'dip' in h'(j2{i') at 2 A, 
particularly in the clean limit. Otherwise the curves all approach the Drude limit at high 
frequency, which, for this property, is unity (conductivities are in units of ne^/m = Up/An). 



6.2.2 Penetration Depth 

Before we examine the effects of the electron phonon interaction on the real and imaginary 
parts of the conductivity, we first summarize the 'BCS' results for the penetration depth as a 
function of impurity scattering, which can be extracted analytically ||186| , |274| from the zero 
frequency limit of the conductivity. The result is, with a = 



X\T = 0) Xl{T = 0)[2a aVT^ 



sin fvl - a'^)> a <1 




Here, the zero temperature London penetration depth in the clean limit is given by 



A^(^-0)^£l. (143) 



In the weak scattering limit Eq. ( |142| ) reduces to the more familiar form 



1 1 1 



A2(0) A2(0)l + fa- 



(144) 



This expression can be written in terms of the zero temperature coherence length, ^q, and 
the mean free path, i, using A = ^ and vp = i/r, where vp is the Fermi velocity: 

111 , , 

(145) 



A2(0) Xl{0)l + ff 
6.2.3 Microwave Regime: Coherence Factors 

The microwave regime (1 - 60 GHz) corresponds to very low energies (1 GHz = 0.0041 
meV). This energy scale is much lower than that of the superconducting energy gap. Mea- 
surements of the microwave response of a superconductor have been used in recent years to 
determine the penetration depth and optical conductivity in the high Tc cuprates ||275| , |276|| , 



but, historically, either the real or the imaginary component of the surface impedance was 
measured, making a determination of the complex conductivity impossible. It is of interest 
to examine the conductivity in this case, because BCS theory makes a highly non-trivial 
prediction that the real part of the conductivity shows a so-called coherence peak just below 
Tc. This coherence peak was almost simultaneously predicted and observed |p77| , |278|| in 



measurements of the NMR relaxation rate ||279|] . We will briefly discuss the source of these 



coherence factors, and return to a description of the microwave conductivity, since a detailed 
discussion of the NMR relaxation rate [|281|| is outside the scope of this review, and the final 



expression relevant to superconductors is a special case of the microwave conductivity. 
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Within BCS theory the transition probabihties between an initial and final state that 
enter the expression for various linear response functions are of the form |]5|,|46[ 



-Pkk' = (^ikMk' =F ^^k^^k')^ (146) 
with Mk, f k the amplitudes that relate quasiparticle operators to electron operators 

and Ek is the usual quasiparticle energy: 



i?k = V(ek-/i)2 + A^ (149) 

where ek is the electron band energy, /i is the chemical potential and Ak is the gap function. 
In Eq. (|146|) the upper (lower) sign corresponds to case I (case II) observables. These 
signs have important consequences for the response, particularly just below Tc. A case in 
point is the electromagnetic absorption; the temperature-dependent result (derived from 
Eqs. (|13|,|133)) in the dirty limit (1/r » A) is ||53| 



+ - '^^^l Cl (E^ - A^y^iil + ty - Ag) - ^^(^ + 

where Aq = A(T) is the temperature-dependent gap function. The second plus sign in 
E[E + u) + Aq which appears in this expression is due to the fact that the electromagnetic 
absorption is a case II observable. In a case I observable this would be a minus sign; it is 
then readily seen that whereas Eq. (150) contains a divergence as z/ — >■ 0, the corresponding 



case I observable would not, as the numerator (coming from the coherence factor given in 
Eq. ( |146D ) would then cancel the density of states factors, which are explicit in Eq. ( |15CI| ), 



and which contain square-root divergences. In both cases the 'freezing out' of excitations as 
the temperature is reduced leads to a low temperature suppression of the response function 
— this is simply a consequence of the gap. On the other hand, near an enhancement is 
expected for type II observables, while, for type I observables, the response is immediately 
suppressed as the temperature is lowered below the superconducting transition temperature. 
In the limit that the frequency is zero, one obtains from Eq. ( |150|) , 

S = ^^l^<-i) - (VrO./(i/ro„, (151) 

which is formally divergent (at all temperatures). The divergence is in fact eliminated in 
practice by anisotropy in the gap or retardation effects. As noted by the second equality, 
this is the expression for the superconducting to normal ratio of the NMR relaxation rate. 
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For a type I observable (like the ultrasonic attenuation) the numerator in Eq. ( p.51| ) has a 
minus sign, so that numerator and denominator cancel, and the remaining integral is trivial. 
One obtains 

= 2/(A(T)), (152) 

where as{n) is the ultrasonic attenuation in the superconducting (normal) state, and / is the 
Fermi function. This is a monotonically decreasing function as the temperature decreases 
from Tc to zero. 



6.2.4 Far-Infrared Regime — Arbitrary Impurity Scattering 

The expressions for the optical conductivity provided in the last three subsections apply only 
in the dirty limit. As already mentioned earlier, a comprehensive expression (for all values of 
elastic impurity scattering), along with a very efficient FORTRAN program, was provided in 
Ref . |p73|| . For completeness, we illustrate here the temperature dependence for two extreme 
cases, close to the clean limit (1/r = 1 meV), and the dirty limit {1/t —>■ oo), in Fig. 29 and 
Fig. 30, respectively. As noted earlier, the optical gap (= 2A(T)) is clearly evident in both 
the real and imaginary part of the conductivity. The evolution from the normal state to the 
superconducting state is clearly evident as well; note, in particular, that in the real part of 
the conductivity, the missing area is taken up as a delta function at the origin (not shown). 



6.3 Eliashberg Results 

Within Eliashberg theory, changes occur for two related reasons. First, even in the normal 
state the self-energy acquires a frequency dependence (no wavevector dependence, because 
of the simplifying assumptions made at the start); secondly, the gap function in the super- 
conducting state acquires a frequency dependence and acquires an imaginary part. This 
latter fact tends to smear many of the 'sharp' results shown in the last section, a feature 
which is already evident in comparing the single electron densities of states in Fig. (24) to 
those in Fig. (21b), for example. For this reason, it is important to re-examine the impact 
of retardation on a variety of observables. 



6.3.1 NMR Relaxation Rate 

In the first few years following the discovery of the high temperature superconductors P], 
several anomalous features were measured in the superconducting state. One of these was 
the absence of the coherence peak (the so-called 'Hebel-Slichter' peak) in the NMR spin 
relaxation rate, 1/Ti, just below Tc ||282| |. Motivated by the possibility that this 'anomaly' 
could be explained by damping effects due to retardation, Allen and Rainer ||168|| and Akis 



and Carbotte [ ^67|| calculated the ratio of the relaxation rate in the superconducting state 
to that in the normal state with several hypothetical electron-phonon spectra (obtained 
by scaling known spectra from conventional superconductors). Both groups found that 
sufficiently strong coupling (as measured by A or Tc/uJ\n) smears out the coherence peak 
entirely. An example is shown in Fig. 31 (taken from Ref. [ |168|| ), which shows the theoretical 
and experimental ||283|| results for a conventional superconductor (Indium) along with data 
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from YBCO |P84|| , and theoretical results obtained using scaled spectra. While the present 
consensus is that the lack of a coherence peak is not solely due to damping effects, the lesson 
learned from these calculations is clear: retardation effects damp out the coherence peak in 
the NMR relaxation rate. It is worth noting here that even within a BCS framework (i.e. no 
retardation), the coherence peak can be suppressed in the dilute electron density limit ||285|| . 

6.3.2 Microwave Conductivity 

A natural extension of this argument applies to the microwave conductivity. In this case, 
even within BCS theory, a divergence does not occur since the experiment is conducted at 
some definite non-zero microwave frequency (see Eq. ( |150| )). Before discussing retardation 
effects, however, it is important to realize the amount of impurity scattering (as characterized 
by 1/r) also influences the height and presence of the coherence peak ||260| , |261| . In Fig. (32a) 



we show, within the BCS framework, the conductivity ratio for a small but finite frequency 
as a function of reduced temperature, for a variety of elastic scattering rates, ranging from 
the dirty limit to the clean limit. Quite clearly the coherence peak is reduced and then 
eliminated function of 1/r. 

To see how retardation effects also serve to reduce and eliminate the coherence peak 
(just as in NMR) we focus on the dirty limit (1/r ^ oo) where the peak is largest without 
retardation. In Fig. (32b) we show results obtained from a Pb spectrum (Fig. 11), scaled 
by varying degrees to increase A from 0.77 to 3.1. For the largest coupling considered the 
coherence peak has essentially vanished. This is the same effect seen in the NMR relaxation 
rate. In Fig. (32c) we illustrate the impact of changing the microwave frequency. Clearly, in 
the limit of very weak coupling (BCS) one expects the strongest variation, since, as z/ — 0, 
the BCS result will diverge logarithmicly. However, as the coupling strength increases, 
the damping due to retardation reduces the peak far more effectively than an increase in 
microwave frequency would, so that the conductivity ratio (at some temperature near where 
a maximum would occur in the BCS limit) is essentially constant as a function of frequency. 
This is clearly illustrated by the two lowest curves in the Figure, representing the strongest 
coupling situations. 

A measurement of the coherence peak in the microwave wasn't actually performed until 
the early 1990 's, in Pb ||286|| and in YBCO ||287|| (although the peaks observed in these latter 



measurements are now thought not to be the BCS coherence peak ||276| , p88| , p8£ 



Several other groups have since examined the microwave response in conventional super- 
conductors. In Ref. p90|| Nb was examined in detail. The experiment was performed at 



17 GHz, and a prominent coherence peak was observed, as shown in Fig. 33. Also shown 
are theoretical curves obtained from Eliashberg calculations; they all fall significantly below 
the experimental results. We have also included the BCS result (dotted curve) computed 
for this frequency; it is not very different from one of the curves obtained using the full 
Eliashberg formalism. The BCS result represents probably the highest achievable coherence 
peak; other alterations of the standard theory (anisotropy, finite bands, non-dirty limit, etc.) 
would tend to decrease the theoretical result further. Hence, at present the coherence peak 
observed in Nb remains anomalous because it is too big. Other measurements in Nb and 



Pb [p91|| showed agreement with Eliashberg theory, but they were carried out at a much 
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higher frequency (60 GHz). Another measurement of the electrodynamic response (using 
simultaneous measurement of the amphtude and phase of the transmission in Nb thin films) 
P92|| supported our results. A more recent measurement of the coherence peak in NbaSn 
|293|| also finds a large discrepancy with Eliashberg theory — the experimental results show 
a peak which is far too large compared to theory. 



6.3.3 Far-Infrared Regime 



While more recent investigations of the far-infrared (and slightly lower Terahertz) regime 
in superconductors utilize transmission techniques which simultaneously measure amplitude 

the more conventional Fourier-transform spectroscopy ||180|| 



and phase information ||288| , 
requires Kramers-Kronig relations, as outlined in Section (3.3.3). For this reason the entire 
spectrum needs to be measured, often with an assortment of spectrometers [p.79|| . How do the 
real and imaginary parts of the conductivity change as a function of the coupling strength A ? 
In Fig. 34 we show real ((a) and (b)) and imaginary ((c) and (d)) parts of the conductivity 
with 1/r = 2 meV and 25 meV, respectively. In all four figures it is clear that an increased 
coupling strength decreases the real and imaginary parts of the conductivity, at least in the 
low frequency regime. In fact, at low temperatures, in the normal state, one can derive a 



Drude-like expression \\IS4 



ne ±11 

where m*/m = 1 + Aand r/r* = 1/(1 + A). This expression clearly indicates that, while the 
zero frequency conductivity remains unaffected, the rest of the conductivity is diminished 
by the electron-phonon interaction [|100|| . In fact integration of Eq. (|153|) yields the result 



1/7 



\ dv (Td 



Tine 



2 m* 1 + A 



(154) 



This is lower than the Kubo sum rule [p57|| by the factor of 1/(1 + A), which says that the rest 
of the area is taken up in the phonon-assisted absorption, which occurs at higher frequency 
(in the phonon range). Also note that one effect of an increased electron-phonon interaction 
strength is to decrease the impurity scattering rate: 1/r — > 1/tj^. This occurs because the 
inelastic scattering reduces the spectral weight of the quasiparticle undergoing the elastic 
scattering. Further discussion of the Drude-like behaviour at low frequency but for non-zero 
temperature can be found in Ref. |p.83| , |184j . 



Returning to Fig. 34, we note that except for small corrections to the gap edge as A 
increases (2 A tends to increase as well), the occurrence of an abrupt onset of absorption in 
the real part ((a) and (b)) exists for all coupling strengths. While a cusp remains in the 
imaginary part ((c) and (d)), its size is clearly diminished as the coupling strength increases. 
Note that the penetration depth (given by the square-root of the inverse of the intercept in 
the imaginary part — see Eq. ( |141[ )) tends to increase as the coupling strength increases. 
Also note that, while not apparent on the frequency scale shown in (c) and (d), the frequency 
times the imaginary part of the conductivity approaches unity (in units of ne^/m) as the 
frequency approaches large values. This fact was utilized in the case of Bao.6Ko.4Bi03, which 
we briefly discuss next. 
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Fig. 35 shows the imaginary part of the conductivity obtained from reflectance measure- 
ments on Bao.6Ko.4Bi03 



186,294,295 



A prominent dip occurs near 12 meV, which has 
been roughly fit by two models as indicated. The occurrence of this dip fully supports the 
existence of a superconducting state with s-wave symmetry, with a gap value that is high 
compared to that expected from BCS theory (2A//cbTc ~ 5 compared with 3.5). This value 



is somewhat higher than that obtained previously with infrared ||296|] or tunneling | 297 , 298 



measurements. Nonetheless, a thorough analysis of the temperature dependence of the Drude 
fits at low frequency [ p.84|| and the frequency dependence illustrated in Fig. 35 ||186|| shows 
that the electron-phonon interaction must be weak in this material, too weak to support 30 
K superconductivity. Two model calculations are shown with the data in Fig. 35. The data 
is clearly consistent with an electron-phonon coupling strength A ~ 0.2 (which requires an 
additional mechanism to produce Tc = 30 K), and entirely inconsistent with A ~ 1. 

As is clear from the preceding paragraph, either the real or the imaginary part of the 
conductivity contains all the relevant information about the absorption processes in the 
system. This is due to the fact that they obey Kramers- Kronig relations, which ultimately 
can be traced to requirements of causality and analyticity | ]18U|| . In an effort to make these 
absorption processes more explicit, one can also favour other functions; a particular example 
is the effective dynamical mass, m*(z/), and the effective scattering rate, l/r(z/), introduced 
through [ p.83|| 

^(-) = #TT^^^— rrr^, (155) 



Air 1/t(z/) — ium* (u) /m' 

where up and m are the bare electron plasma frequency and mass, respectively. Then, 
one can define an effective scattering function, l/r(z/), which can be extracted (say, from 
experiment) through 

' (156) 



l/r(.) ^ ^Re4^ 



bearing in mind that cr(i/) itself has been obtained through Kramers-Kronig relations from, 
say, reflectance data. This is precisely the function required to invert normal state conduc- 
tivity data to extract a^F(z/) (see Eq. (^)). A plot of 1/t(z/) vs. frequency is neverthe- 
less revealing. It tends to illustrate at roughly what energies absorption process 'turn on' 
183| , |299|| . For example, we show in Fig. 36 the function l/r(i/) derived from conductivity 
results of model calculations for Bai-ajK^BiOsand YBaCusOr-^; [|181|]. The former uses a 



model phonon spectrum extracted from neutron scattering measurements [pOQI while the 
latter uses a model spin fluctuation spectrum |p01|| . The fact that the YBaCusOy.a; result 
continues to rise at 300 meV reflects the frequency scale of the spin fluctuation spectrum. In 
contrast, the Bai-ajK^jBiOs result has almost saturated by 100 meV, since the phonon spec- 
trum extends only to 80 meV. More detailed comparisons with self-energy- derived scattering 
rates have been provided in Refs. ||181| , |183| . 



The results shown in Fig. 36 were obtained in the normal state. In the superconducting 
state the presence of a gap will modify the low frequency behaviour of the scattering rate, 
l/r(i/). Results within BCS theory (elastic scattering rate only — no inelastic scattering) 
are shown in Fig. 37. At low frequencies the overall scale of the effective scattering rate is 
set by the elastic scattering rate (2 and 25 meV, respectively). Note that in the gap region 
(below 32 meV) the effective scattering rate is zero (at zero temperature), while slightly 
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above the gap the effective scattering rate below is actually enhanced with respect to the 
normal state value. In Fig. 38 we show the effective scattering rate vs. frequency using 
the model Bai-j^K^-BiOs spectrum in (a) the clean limit and (b) with significant impurity 
scattering. The results are qualitatively similar to those in Fig. 37. 



6.4 Phonon Response 

Much of this review has focused on various properties whose determination allows one to infer 
the degree of electron-phonon coupling that exists in the material under study. The majority 
of properties that fall in this category refer to a modification of the electronic structure or 
response due to a coupling with phonons. To a much lesser extent the phonons themselves 
are modified because of the electron-phonon coupling, and in this section we briefly address 
a few examples in this category. 

The impact of the superconducting state on the phonons was first investigated using 
ultrasound experiments ||302|| . Sound waves are attenuated due to their absorption in the 



solid. The absorption requires interaction with electrons with energies very close to the 
Fermi energy (the phonon energy is typically very low for sound waves — in the 100 MHz 
= 0.0004 meV range). These electron states are gapped in the superconducting state, so the 
attenuation is expected to be suppressed to zero as T — 0. The BCS result, given by Eq. 
( |152|) , is valid for an order parameter with s-wave symmetry. A similar law can be derived 



for other symmetry types [303], which results in some sort of power law decay rather than 
exponential at low temperatures. 

Of main interest here is how Eq. ( |152| ) is modified when retardation effects are ac- 



counted for. An early calculation ||304|| found that retardation effects did not alter the result 
Eq. ( |152| ). Therefore, little can be learned about the electron-phonon interaction through 
ultrasonic experiments; instead, one should examine higher energy phonons. 

The classic experiment of this type was performed using neutron scattering on Nb and 



NbaSn |P05|| . The idea is simply that the electron charge susceptibility modifies the phonon 
spectrum. Within the normal state this modification is hardly noticeable in metals over a 
temperature range of 300 K or so. However, when the material goes superconducting, the 
electron density of states is profoundly modified at energy scales of order the gap; this in turn 
will affect phonons whose energy is on the same scale. In particular, a low energy phonon 
(energy less than 2A) that had a finite lifetime because it could decay into an electron-hole 
pair will be unable to do so in the superconducting state because no states exist at energies 
below the gap, A. Therefore its lifetime will lengthen considerably in the superconducting 
state, resulting in a narrower lineshape below Tc. Fig. 39 shows the experimental result from 
NbsSn ||305|| where the lineshape has clearly become narrower in the superconducting state. 
Similarly, if the phonon energy is slightly above 2A, then, under the right conditions, the 
linewidth will increase, since the electron density of states increases in this energy regime in 
the superconducting state. 



A detailed theory of these effects was first given in Ref. ||306|| , within BCS theory. The 



theory consists of a calculation of a response function corresponding to a Case I observable. 
Similar calculations were performed much later by Zeyher and Zwicknagl |268] to understand 



the frequency shifts and linewidth changes (due to superconductivity) in the q = Raman 
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spectra for various optical modes in YBaCuaOr-^:. They found, using the BCS approxima- 
tion, 



Re An(q = 0,z/ + z,5) _ -jTHF^^^ (t!^) for z7 < 1 

iV(0) [ —^in(2u'^ - 1 + 2vVW^) for u > I. 

The imaginary part is given for all temperatures by: 



(157) 



Jm An(q = 0, z/ + 25) tanhpp/A 

where v = z//(2A(T)). Here, An(q, u + i5) is the change in the phonon self energy between 
the superconducting state and the normal state. A positive (negative) real part means that 
phonons harden (soften) in the superconducting state, while a positive (negative) imaginary 
part means that the phonon linewidths narrow (broaden). Thus, phonons below the gap 
edge (2A) soften while those above harden. Also, above the gap edge they broaden while 
below their linewidth does not change. The broadening above 2 A can be understood as being 
due to the enhanced scattering with electrons, since the electron density of states now has a 
square- root singularity in the energy range of A, and the phonon self energy is essentially a 
convolution of two single electron Green functions (see Eq. ( p.28| )). 

Eqs. ( |157|J158|) have been derived assuming single particle Green functions without im- 



purity scattering. The q = limit is somewhat anomalous in this case, in that the phonon 
width is already zero in the normal state. Hence, no change can occur in the linewidth in the 
superconducting state, for frequencies below 2 A. A calculation with impurities | |269| | provides 
a non-zero linewidth in the normal state. Because of the gap in the single electron density of 
states in the superconducting state, this linewidth is reduced to zero when the system enters 
the superconducting state, so the change in the imaginary part of the phonon self energy is 
positive. These results are summarized in Fig. 40. Note that the softening below the gap 
edge is significantly reduced with impurity scattering present, and the phonons above 2 A 
also soften when a significant degree of impurity scattering is present. As Fig. 40b shows, 
phonons whose energy lies below 2A acquire a narrower linewidth in the superconducting 
state, as noted above. 

The effects of retardation on the phonon self energy are not very significant. The changes 
that do occur follow the changes already discussed due to including elastic scattering; high 
energy phonons soften rather than harden, and the broadening that accompanies this soft- 
ening is reduced compared to the clean BCS case. More detailed changes are documented 
inRefs. |26|,|69l. 



Because these phonon changes can be observed through neutron scattering experiments. 



it is of interest to examine the phonon self energy at non-zero momentum, q ||270| , |27l[| . In 
this case the phonon has a non-zero linewidth in the normal state, and so line narrowing 
is observed in superconducting state at low frequencies, due to the development of a sin- 
gle electron gap. The detailed frequency dependence is a function of the band structure; 
in particular, with two dimensional nesting phonon changes due to superconductivity are 
enhanced [p71|] . 
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7 Summary 



We have examined a variety of ways in which the retarded electron phonon interaction 
influences the properties of a conventional superconductor. The first and simplest effect is 
through a renormalization of Fermi Liquid parameters, like the effective mass. While this 
effect appears in a number of normal state properties (for example, the low temperature 
electronic specific heat capacity, where the Sommerfeld 7 is enhanced by 1 + A — see Eq. 
( p.2CI| )), it also appears in many superconducting properties. The most obvious (but least 
measurable) example is in the Tc equation, Eq. (p2D, where 1 + A appears in the exponent. 
Another (perhaps more detectable) occurrence is in the slope of the upper critical magnetic 
field. In each of these cases, the renormalization occurs in the normal state — its occurence 
in the superconducting state is because the property in question depends on the normal state 
effective mass, or Fermi velocity, etc. One should also bare in mind that the factor 1 + A, 
comes from a weak coupling approach. In a strong coupling approach, an electron phonon 
renormalization is still present, but may be much more significant than suggested by the 
weak coupling approach, and polaron-like physics may dominate []97 |. 

The most important manifestation of the electron phonon interaction is the supercon- 
ducting state itself. In fact, according to our present understanding of Cooper pairing, the 
electron phonon-induced attraction between two electrons would not overcome their direct 
Coulomb repulsion, except for the fact that the former is retarded whereas the latter is not. 
This gives rise to the pseudopotential effect; in some sense the pseudopotential effect is the 
true mechanism of superconductivity, rather than the electron phonon interaction per se. 
This is perhaps emphasized in the cuprate materials, where presumably the electrons could 
not utilize the difference in energy (and hence time) scales between the attractive mecha- 
nism (whatever it is) and the direct Coulomb repulsion to overcome the latter. Instead the 
pairing has apparently adopted a different symmetry (d-wave) to avoid the direct Coulomb 
repulsion. 

Nonetheless a minimal accounting for these retardation effects accounts fairly well for 
the superconducting ground state. This was accomplished by BCS theory. A more accurate 
theory with retardation effects (Eliashberg theory) quite clearly accounts for quantitative 
discrepancies with experiment. Here, Pb and Hg are held up as paradigms for retardation 
effects, the simplest occurring in a measurement of the gap ratio, for example. The BCS 
theory predicts a universal number for this ratio, 2A/ kBT^ = 3.53. With Eliashberg theory a 
value for Pb is found close to 4.5, in excellent agreement with experiment. We have charac- 
terized the discrepancy with BCS theory through a retardation parameter, Tc/ui^. Various 
properties have been quantitatively accounted for through simple analytical expressions with 
this parameter, as given in Sections 4 and 5 (see Ref. |Tl|] and references therein for many 
more) . 

Finally, various dynamical properties exhibit 'signatures' of the electron-phonon pairing. 
These tend to manifest themselves as 'wiggles' in the data, the most famous of which occurs 
in the tunneling data, and allows an inversion to extract the electron phonon spectral func- 
tion, a^F(z/). As we saw briefly in Section 3, and then again in Section 6, these 'wiggles' 
occur in various two-electron response functions, most prominent of which is the optical 
conductivity. An accurate measurement of these response functions allows one to infer a 



56 



significant electron-phonon coupling. 

We have focussed on very conventional superconductors, and have, for example, avoided 
any analysis of the high temperature superconductors. Signs of electron phonon interactions 
have occurred in these new materials as well, but the relation to the superconductivity in 
them is yet unclear. Moreover, such effects will no doubt be covered in other chapters. 
Nonetheless, we wish to add a few remarks about other classes of superconducting materials 
that have been discovered over the last twenty years. 

Cubic Perovskites, beginning with strontium titanate (SrTiOs) [^1],^, have already been 
discussed in Section 2. As mentioned there, these compounds (including BaPbo.75Bio.25O3 
(Tc ~ 12 K) and Bai.^K^BiOa (T^ ~ 30 K) |^) are generally regarded as in a distinct 
class from the high cuprates. This has left them, somewhat by default, as electron-phonon 
driven superconductors. On the other hand, there is strong optical evidence ||186| , |294|| that 



the electron phonon interaction is very weak in these materials. Hence, as far as we are 
concerned, the mechanism of superconductivity in these perovskites is not understood at all. 
Tunnehng studies ||297| , |298|| are divided on this issue. 

One- and two-dimensional organic superconductors were discovered in 1979 [p07|]. The 



subject had developed sufficiently so that, by 1990, a book devoted to the topic was written 
| |308| |. Organic superconductivity represents another interesting idea that was first presented 
by theorists [|106| , |107|| , on the basis of a phonon-mediated interaction, but that now is con- 
sidered by most practitioners not to be due to electron phonon interactions. Many of the 
organics abound in physical phenomena, with several containing, on the same phase dia- 
gram, charge density wave (CDW) and spin density wave (SDW) instabilities, juxtaposed 
with superconductivity |P09|| . The nature of the superconducting state has not really been 
sharply defined by experiments, to the extent that both singlet and triplet pairing may be 



present |piO|| , and the presence of a gap has not been unequivocally established. While it 
is probably fair to say that the electron phonon interaction has not been ruled out as the 
mechanism for superconductivity, spin fiuctuation-mediated pairing seems to be favoured 

lua. 



Heavy Fermion systems were discovered to be superconducting also in 1979 [311]. While 



Tc has remained low, these compounds have remained of interest because (i) the root cause of 
the heavy electron mass is not completely understood, and (ii) the superconducting ground 
state coexists in a number of cases with antiferromagnetic order. It has now been established 



through thermal conductivity measurements that the order parameter contains nodes [312 



and the circumstantial evidence points towards an unconventional magnetically mediated 
mechanism for superconductivity [pl3|[ . There is very little indication that superconductivity 
in this class of compounds has anything to do with the electron-phonon interaction. 

Superconductivity in alkali-doped buckminster fuUerene (A3C60, with A = K, Rb, Cs) 
was briefiy mentioned earlier in this chapter. On the basis of optical measurements ||314|| , 
a sizable electron phonon coupling was inferred, and, in fact q;^F(z/) was extracted by an 
inversion procedure outlined in Section (3.3.3) [|178|[ . Evidence for electron phonon-mediated 



superconductivity was also presented in earlier reviews [|7^. On the other hand, doubts 



remain concerning the validity of a weak coupling framework |pl5|[ . One would like to 



understand the 'bigger picture', i.e. the progression from insulator with pure Cgo through 
the superconducting phase with A3C60, and back to insulator with AgCgo- In fact, band 
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structure calculations ||316|| suggest (simplisticly) that A2C60 should be superconducting 
with a higher than A3C60, when, in fact, that compound does not readily form. 

The electron phonon theory can be subjected to even more tests, now that workers have 
managed to fabricate a field effect transistor which allows electron ||317|| and hole ||318|] doping 
of Ceo- Tc is much higher for hole doping (Tc = 52 K), and spans a wide range of dopant 
concentration. In fact this peculiar asymmetry between electron and hole doping finds a 

rmi 



natural explanation through the hole mechanism of superconductivity | pi9| j. An explana- 
tion in terms of a dopant-dependent electron phonon coupling strength appears somewhat 
unnatural. 

The borocarbides (RNi2B2C, where R denotes a rare earth element) were found to be 
superconducting in 1993 ||320 , 321 1. In addition to having a sizeable transition temperature 
(Tc «i 20 K), some of these compounds exhibit coexistent superconductivity and antiferro- 
magnetic order, and indeed, share some similarities with the heavy fermion compounds |p22 



Nonetheless, tunneling has determined that a well-defined gap exists at low temperatures, 
and this and other measurements have established these compounds to have very BCS-like 
properties | |323| j . A detailed comparison of various superconducting properties with results 
based on Eliashberg theory (including some small anisotropy) ||324|| yields excellent agree- 
ment. A model spectrum was used for the electron phonon interaction, and, at present, it 
remains unclear to what extent this agreement points unequivocally to the electron phonon 
mechanism for superconductivity in these compounds. 

Very recently, superconductivity with Tc = 39 K has been discovered in the very simple 
binary compound, MgB2 ||325|1 . Preliminary results indicate a gap in the single electron 

Calculations of the 



density of states [ 326 - p28|l , and an isotope effect has been observed [ 32£ 
electron phonon coupling strength, not quite consistent with Tc = 39 K, have been reported 
| 330|| , as has a competing non-electron phonon mechanism, based on the hole mechanism 
||331|| . More experimental results will be required before a real assessment of the electron 
phonon mechanism can be provided. 

Finally, Sulfur has been found to exhibit a high superconducting transition temperature 
(Tc = 17 K) ||332|| . Very little work has been carried out regarding the mechanism; a notable 
P33l. 



exception is Ref. |P33|| , where ab initio calculations are performed to estimate the electron 
phonon coupling strength for Sulfur. They find that under pressure, in a different structural 
phase, the electron phonon coupling is enhanced, consistent with the increase in Tc. 

As is evident by the foregoing examples, a steady search for new superconductors is 
being rewarded with discoveries of materials with high critical temperatures, now in the 
same category as those of the high temperature cuprates. The A15 compound record of 
Tc ~ 23 K would have been broken many times by now, even if the layered cuprates had 
not been discovered. Most intriguing is the fact that many of these compounds may be 
driven to the superconducting state through the electron phonon mechanism. As far as 
future developments in this area is concerned, an obvious question to be addressed is the 
soundness of the original Cohen- Anderson estimate [|194|| for the maximum electron phonon 
mediated critical temperature. It may simply be a matter of quantitative assessment, or 
perhaps some more exotic effect (within the electron phonon picture) has been overlooked. 
An intermediate or strong coupling approach [^] may yet provide new insights. Finally, one 
can't help but notice the recent resurgence of investigations in the high temperature cuprates 
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themselves, that indicate strong electron phonon effects | 334 |. To paraphrase | 335 |, 'The fat 
lady probably hasn't yet sung'. 



8 Appendix: Microscopic Developments 

In this Appendix, we will first outline a derivation of Eliashberg theory, based on a weak 
coupling approach. By this we mean that we start with momentum eigenstates. While other 
derivations may be given in other chapters, we include one here to keep this chapter some- 
what self-contained. Migdal theory follows by simply dropping the anomalous amplitudes 
in what follows. We will then outline various other attempts to understand electron phonon 
interactions, particularly in the strong coupling regime. 



8.1 Migdal-Eliashberg Theory 

We begin with the definition of the one electron Green function, defined in momentum space, 
as a function of imaginary time |^ , 

^(k, r-T') = -< T,Ck.(r)cL(r') >, (159) 

where k is the momentum and a is the spin. The angular brackets denote, as usual, a 
thermodynamic average. With this definition such a Green function can be Fourier expanded 
in imaginary frequency: 

-| oo 

G(k,r) = -^e-^'^™-G(k,z^^) 

P -oo 

G{\i,ium) = f dTG%T)e"^^\ (160) 
Jo 

The frequencies ium are known as the Matsubara frequencies, and are given by ium = 
27rT(2m — 1), m = 0, ±1, ±2, where T is the temperature. Because the c's are Fermion 
operators, the Matsubara frequencies are odd multiples of ivrT. As is evident from these 
equations, the imaginary time r takes on values from to /? (= -j^)- 
Similar definitions hold for the phonon Green function: 

D(q, r - r') = - < T,A^{t)A_^{t') >, (161) 

where y4q(r) = aq(r) + at.q(r). The Fourier transform is similar to that given in Eq. ( |160D 
except that the Matsubara frequencies are ivn = i7TT2n, n = 0, ±1, ±2, ... i.e. they occur at 
even multiples of inT. 

To derive the Eliashberg equations, we use the equation-of-motion method, taken from 
Ref. The starting point is the (imaginary) time derivative of eq. ( |159| ) 

^G(k,r) = -(5(r) - <r,[/J-/iAr,Ck.(r)]cL(0)>, (162) 
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where, without loss of generahty, we have put r' = 0. For definiteness, we use the Hamilto- 
nian (|^), and, in addition, assume, for the Coulomb interaction, the simple Hubbard model, 
Hcoui = UY,ifT-i^fT-ii- The sum result is 

kcr q VjV kk' k,k',q 

(163) 

where the various symbols have already been defined in the text. Working out the commu- 
tator on Eq. ( |162| ) is then straightforward. We obtain 

+ ek^G't(k,r) = -5{t) - -^^^kk' < 7;v4k„k'Mck'TW4T(0) > 



U 

pp' 



+ ]V E < ^rC;.,_k+pi(^)cp'i(^)cpT(^)4t(0) >, (164) 



where for definiteness we are considering the Green function with a =]. On the right- 
hand side of Eq. ( |164| ) various higher order propagators appear; to determine them an 
equation of motion would have to be written, which would, in turn, generate even higher 
order propagators, eventually leading to a set of equations with hierarchical structure. This 
infinite series is normally truncated at some point by the process of decoupling, which is 
simply an approximation procedure. For example, in Eq. ( |164| ) the Coulomb term is normally 
not expanded further; instead a decoupling procedure is employed. Thus, under normal 
circumstances, the last term would become 

< ^rCj,,_k+pi(^)Vi(^)CpT(^)CkT(0) > ^ < ^rCi,,_k+pi(^)Cp'i(^) >< ^rCpt (r)4^ (0) >, 

^ -(5kpGi(p',0)GT(k,r). (165) 

The case of the electron-phonon term is a little more subtle, however. In this case we define 
a Green function, 

G2(k,k',r,ri) T.Ak„k'(^)ck'T(n)4^(0) >, (166) 
and write out an equation of motion for it. We get 



^G2(k, k', r, n) = -^k-k' < T,Pk-k'(r)ck'T(ri)4^(0) >, (167) 



where Pqij) = aq(T) — a_q(r). Taking another derivative yields 



-CJk-k' G'2(k,k',r,ri) = XI 2t^k-k'fi'k-k' < ^r4//_k+k'a(^)Ck"<7(^)Ck'T(n)4^(0) > . 

k"cr 

(168) 

One might be tempted to decouple Eq. (|168D and thus close the hierarchy that begins 
with Eq. (|164|) . However, retardation effects are properly included only when the phonon 
propagator is taken into account. While the electron-phonon interaction affects the phonons 
as well as the electrons, the influence on the phonons occurs most at higher temperatures. For 
many materials the phonons have reached their ground state configurations by about room 
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temperature. As a result, for low temperatures the phonons remain virtually unaffected by 
the electron-phonon interaction, and it suffices to disregard the electron-phonon interaction 
as far as the phonons are concerned provided they have been properly renormalized due to 
effects which took place at higher temperature. To put this another way, inelastic neutron 
scattering measurements of the phonon dispersion curves show a dependence on temperature 



only at temperatures well above room temperature ||89| , 

As already mentioned in the text, the phonons are normally taken from experiment, and 
hence the "calculation" of the phonon propagator is greatly simplified. One simply assumes 
that the phonons are non-interacting. The equation of motion for the phonon propagator is 
then 




\D{q,T-T') = 2ujq6{T-T'). (169) 
Utilizing this expression in eq. (|168|) then yields 



G2(k, k', r,T) = ^Y.f dT'g^^,D{k - k', r - r') < T,4„_k+k'<x(^')ck'v(r')ck'T(r)cLt(0) >, 

(170) 

where now ti has been set equal to r as is required in Eq. (|164|) . This can now be substituted 
into Eq. ( |164| ), and the whole result can be Fourier transformed (from imaginary time to 
imaginary frequency). Before stating the result of this exercise, however, we note that the 
superconducting state is specially characterized by the existence of anomalous amplitudes, 
attributed to Gorkov [Q] and often referred to as Gorkov amplitudes. Thus, in the Wick 
decomposition of the various two-particle Green functions, the anomalous amplitudes 



also must be taken into account, in addition to the normal amplitudes given, for example, 
in Eq. (|l65|). 

The anomalous amplitudes take the form 

F(k,r) = - < T.CkT(r)c_ki(0) > (171) 

and 

F{k,T)^-<Tj_^^{T)cliO)>. (172) 

Now it is necessary to go through the same procedure with F and F as with G. The 
methodology is the same, so we skip the necessary steps. 

We then define two self-energies, the usual one (generalized to the superconducting 
state), denoted by S(k, ium), and an anomalous self-energy, often called the pairing function, 
(j)(k,iujrn), and we arrive at Eqs.(H3-H0). 



8.2 The Polaron Problem 

A rather different and less developed approach to the electron phonon problem focuses on 
the effect of the phonons on a single electron. A review is provided in Ref. [0, and we 
merely highlight some of the important points here. 

There are many kinds of polarons, i.e. small vs. large, weakly coupled vs. strongly coupled, 
Frohlich vs. Holstein, etc. As far as we can tell these classifications are merely qualitative. 
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so that, in most cases, distinctions can be readily drawn for extreme parameters only. A 
case in point is the distinction between an itinerant vs. self-trapped polaron. It seems clear 
that no such transition exists, but nonetheless a crossover occurs to a regime in which the 
polaron acquires a very large effective mass. 

In thinking about the polaron problem, there is the usual competition between kinetic 
energy (measured by the hopping integral, t, or the bandwidth, D = 2zt, where z is the 
coordination number for a cubic lattice [z = 2,4 and 6 in 1,2, and 3 dimensions, respectively)) 
and the potential energy (measured hy g — see Hamiltonian (|^)). In addition the phonon 
frequency represents a third energy scale. In the case of the Holstein model, Eq. (^, this 
scale is conveniently represented by a single number, the Einstein oscillator frequency, ue- A 
dimensionless coupling constant, A = 2g'^ / {Duoe)-, corresponds roughly to the enhancement 
parameter introduced in section 2.3 (see Eq. (pTD). Note that in terms of the parameters 
of the original Holstein Hamiltonian, Eq. (|^), A = a^/{KD). An increase in A signifies an 
approach to the strong coupling limit. On the other hand the adiabatic (anti-adiabatic) limit 
is represented by ujE/t —>■ 0{oo). The values of both ratios strongly influence the number of 
phonons present. An early review that clearly delineates these different regimes is provided 
byRef. |33§. 

There have been many approaches to solving the polaron problem (as governed by a 
Hamiltonian like Eq. (^ ) . Some of the early techniques are amply covered in Ref . ||5^ ; these 
are exemplified by weak and strong coupling perturbation theory, and variational methods. 
A review of the perturbation approaches is given in Appendix B and C of Ref. ||337|| . Weak 
coupling follows the Migdal approach, while strong coupling utilizes the celebrated Lang- 
Firsov transformation. 

This transformation immediately results in a narrow band, with effective hopping pa- 
rameter, tes = texp (— (yf^/a;^), along with exponential increases in effective mass and, in the 
adiabatic regime, number of phonons in the ground state ||336| , |338 |. 



With the advent of considerable computing capabilities over the last two decades, exact 
methods have been used, that, in various cases, can span the entire parameter regime. 
The first is Monte Carlo for a single electron, pioneered in Ref. ||339||. Trugman et al. 



161| , |340| , |341|] utilized exact diagonalizations based on a variational Hilbert space obtained 



from repeated applications of the Hamiltonian on a trial state vector; their most recent results 
are capable of achieving very high precision. In the meantime, Proetto and Falicov |p42|| and 
Ranninger and Thibblin [343] used a truncated Hilbert space for a two-site problem, and 
performed a straightforward numerical diagonalization. This was followed by work on larger 
(one-dimensional) lattices (for one electron) in Refs. ||337| , |344| , |345|| . Most of this work was 
performed for a specific model — the Holstein model of electron phonon coupling, already 
referred to in the text. Further work was carried out also for the BLF (SSH) model, in Ref. 
346 1 . Yet another technique utilizes the density-matrix renormalization group (DMRG) ||348|| 
method |p49|| , which has also been extended to many electrons ||350|| . Another variational 



technique known as the Global-Local variational method ||347|1 also provides very accurate 
results for the polaron problem. Finally, two new Monte Carlo methods ||352| , p53|| appear to 
be particularly powerful in obtaining polaron properties. 

In all clearer understanding is emerging; there is no self-trapping transition, in 

any dimension, although there is a farily abrupt (but still smooth) crossover from weak 
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coupling- like to strong coupling-like. This crossover has now been investigated in 1,2, and 3 
dimensions as well as with dynamical mean-field theory, which is exact in infinite dimensions 

An actual transition can be observed in higher dimensions in the adiabatic limit {uje = 0) 
However, this hmit is regarded as somewhat pathological, and not representative of 



354|. 



the general case ||347|| . 

Finally, some work has been performed on the bipolaron problem, i.e 

m 



whether two 
Various discussions of 



polarons bind or not. Much of this work is summarized in Ref. 
the supporting evidence and difficulties of these theories can be found in Refs. [ ^55| - |357| . A 
related problem has been asked and partially answered in Refs. |p.61j , P37| , |358|| , which is: to 
what extent will two electrons interacting through phonon exchange and Coulomb repulsion 
form a Cooper pair, particularly as the stength of the Coulomb repulsion is increased well 
in excess of the effective strength of the attractive phonon-induced interaction ? In other 
words, to what degree does the pseudopotential effect play a role in pairing ? 

In Ref. |p37|| one of us found that pairing persists even when the Coulomb interaction 
strength exceeds that of the electron phonon attraction, and in Ref. ||161|| this statement 
was made more precise (see also Ref. [|359|| ). In particular, binding persists only up to a 
point; for sufficiently large Coulomb repulsion, the pair is no longer bound. While more 
work is required, this finding imphes that the usual pseudopotential reduction, given by 
Eq. (^), may be too strong. Eq. (|75|) , for example, achieves a large reduction in the limit 
I^{Ef) — * oc, whereas the result of Ref. [|161|| says that for two electrons, at least, the binding 
is lost in this limit. 



8.3 Many Electrons on a Lattice 



The problem of many interacting electrons is, in many ways, significantly more difficult than 
that of one or two electrons. The dimension of the Hilbert space grows exponentially, so that 
exact diagonalizations become prohibitive. A review of methods and results can be found 
in Ref. [BSOj. As far as the electron phonon problem is concerned, there is some limited 



work which utilizes direct diagonalization, usually in the context of the t-J model |361, 362 



Mainly, however, this problem has been approached through Monte Carlo methods, and a 
variety of (somewhat uncontrolled) Green function techniques. 

Monte Carlo methods have an illustrious history ||363|| . While they are not formally exact 
(because, for example, of a Trotter | |364| | breakup), the error introduced by such a decom- 
position can be controlled. Hence, in principle, and even in practice through extrapolations, 
one can obtain results which are exact to within some known error. 

Some of the first papers to utilize Monte Carlo methods in many body fermion problems 
(in the condensed matter context) addressed the electron phonon problem |365- ^67|] . This 
particular methodology integrated out the fermion degrees of freedom analytically, leaving 
the boson degrees of freedom to which Monte Carlo algorithms were applied. Various mod- 
ifications immediately arose, and were used to address the same electron phonon problem 
| 368| - p7l | ] as well as electron-electron problems ||372|| . Much of this work is reviewed in Ref. 
||68|| ; a more comprehensive review of the many variants of the Monte Carlo method (in 
condensed matter) is provided in Ref. . 



63 



Studies in two dimensions became more feasible in the late 1980's; an immediate question 
that was addressed was the competition between superconductivity and the charge density 
wave (CDW) instability p^, p5| , |373| . At half-filling (where simulations are easiest) the CDW 



instability overwhelms the tendency towards superconductivity, in part because the tight- 
binding model with nearest neighbour hopping exhibits nesting at half-filling. Vekic et al. 
374|| explored the impact of next nearest neighbour hopping (to remove the nesting) but 



found it was difficult to discern whether an incommensurate CDW instability or supercon- 
ductivity dominates. 

Another means of eliminating the CDW is through doping; again most of the work is 
inconclusive. A third means is through the use of a Hubbard U . A study |P75|| of the 



so-called Hubbard-Holstein model in two dimensions found that both the CDW and the 
superconductivity susceptibilities are suppressed as U grows. To our knowledge, however, the 
pseudopotential effect (where the U would essentially cancel the electron phonon interaction 
as far as the CDW was concerned, but not as far as superconductivity was concerned) has 
never been detected in many-electron Monte Carlo studies. 

One of the reasons for exact studies of these lattice models (on small lattices) is for 
use as a benchmark to which diagrammatic methods can be compared. Thus, for example, 
the conclusion in Ref. |^ was that the Migdal formalism, without vertex corrections, 
described the Monte Carlo results fairly accurately, provided phonon renormalization was 
taken into account. A model system is required to determine this, since, in real systems, the 
phonons are often taken from experiment, and already contain renormalization effects. This 
conclusion was confirmed in Ref. ||375|] , as well as in Ref. [|376|| . In this latter reference, the 



authors developed the formalism even further to accommodate a CDW gap, and found good 
agreement with Monte Carlo results. 

Nonetheless, it is probably safe to say that a reliable formalism has not yet been developed 
to investigate low temperature properties of electron phonon systems, particularly slightly 
away from half-filling. A number of attempts have been made, particularly in the case 
of electron-electron interactions [p77| , 378 |, although a comprehensive treatment has not yet 



been achieved (the many-body approaches are also becoming almost as numerically intensive 
as the Monte Carlo methods, and so one of their advantages is diminishing). 

Attempts have also been made to incorporate specific kinds of corrections to the Migdal- 
Eliashberg formalism. One of these categories is the inclusion of vertex corrections. Many feel 
that they may be necessary because the adiabatic ratio uj£,/Ep is not small in some cases 
(eg. high Tc cuprates, and doped buckyballs). In the cuprate materials two-dimensional 
effects may enhance vertex corrections as well. Calculations showing an enhancement of 
due to vertex corrections have been reported in Ref. |]379|| (for a two-dimensional gas). In 



c 



Ref. [380 1, a different tack is taken; Tc is kept fixed, and calculations with vertex corrections 



included can mimic those without through an adjusted /i* (except for the isotope effect). In 



Ref. ||381|| a two-dimensional tight-binding model is used and once again, the conclusion is 
that vertex corrections enhance the pairing interaction. To our knowledge, however, these 
effects have never been observed in exact or controlled calculations. 
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Figure 1: In (a) one electron polarizes the lattice (indicated by dashed circles displaced 
towards uppermost electron); in (b) that electron has moved away. In the meantime a 
second electron (seen below in (a)) is attracted to the polarized region, which has remained 
polarized long after the first electron has left the region. Figure is schematic only, and does 
not, for example, properly convey the opposite momenta such a pair should possess. 
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(b) BLF 
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Figure 2: Schematic of ionic displacements in (a) the Holstein model, and (b) the BLF 
model. In (a) neighbouring chains are distorted in the vicinity of the electron, and in (b) 
neighbouring ions, when displaced while undergoing oscillations, lead to an increased (or 
decreased) overlap region (shaded in black), which leads to an altered hopping amplitude for 
the electron. 
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Figure 3: Real and Imaginary parts of the electron self energy in the normal state, for an 
Einstein spectrum (A = 1). The dotted lines arc the inverse non- interacting electron Green 
functions, a; — (ek — //), for (ek — ^j)/uje — —4, 0, and 4, from top to bottom, respectively. 
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Figure 4: The spectral function for an electron interacting with phonons (Einstein spectrum 
with A = 1) for various momenta as labelled. Note that for each momentum there is a 
delta function contribution (artificially broadened in this figure) whose weight diminishes as 
one moves away from the chemical potential and whose frequency approaches the Einstein 
phonon frequency. The incoherent component grows with increasing ek — jj,, and approaches a 
reasonably well-defined peak centered around ey_ — ii for large values (eg. dot-dashed curve) . 
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Figure 5: Gap anisotropy for Pb as a function of angle, 6', for three different values of 
azimuthal angle, 0. Regions where the Fermi surface of Pb does not exist are indicated by 



vertical lines. Figure reproduced from Ref. ||122|| 
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Figure 6: Normalized tunneling conductance reduced to the corresponding BCS expression, 
for (a) a peak, and (b) a valley in the electron density of states (solid curves). The dashed 
curves were obtained with a constant density of states. Reproduced from Ref. ||134| . 
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Figure 7: Effective electron phonon spectral functions obtained by the inversion of the 
calculated normalized tunneling conductances within the usual (i.e. constant electron density 
of states) Eliashberg theory (solid curves). The input spectral functions are shown with 
the dashed curves. The insets contain the corresponding electron densities of states used. 
Reproduced from Ref. 
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Figure 8: A set of 'constant q' scans in Pb taken at various points along the diagonal in the 
Brillouin zone. Reproduced from Ref. . 
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Figure 10: The electron-phonon spectral function a^F{uj) (solid curve) for PbAoTi.eo deter- 
mined from tunneling experiments and convoluted by instrument resolution of the neutron 
spectrometer compared with the neutron results for the phonon frequency distribution F{uj) 
(dashed curve) measured by incoherent inelastic neutron scattering ||145|| (upper frame). The 
lower frame shows the tunneling results (solid curve) compared with the phonon frequency 
distribution (dashed curve) determined from a Born von Karman analysis of the phonon 



dispersion curves in Ph^oT^.m [|144 
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Figure 11: The electron phonon spectral function a'^F{uj) measured in tunneling experiments 



(dotted curve) compared with that which is calculated from first principles (solid curve) ||148 . 
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Figure 12: I-V characteristic of a Pb-I-Pb junction showing the construction used to find the 
energy gap. The sohd fine and open circles are the current in the normal and superconducting 
states, respectively. Reproduced from Ref. . 



95 





- 




























L 


J 
J 


\ 














1 

■ 




























r 1 
































w 






























r 












w 

J 














— V 





Figure 13: Conductance dl/dV of a Pb-I-Pb junction in the superconducting state nor- 
malized by the conductance in the normal state vs. voltage. Also shown is the two- 
superconductor conductance calculated from the BCS density of states which contains no 
phonon structure. Reproduced from Ref. [52 . 
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Figure 14: Electronic density of states N{E) vs. E — Ao for Pb, obtained from the data 



Fig. 13. The smooth curve is the BCS density of states. Reproduced from Ref. |52[ . 
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Figure 15: The predicted (solid curve) normalized density of states in Pb as a function of 
energy uj compared with measured values (open dots) as a function of energy measured from 
the gap edge. The measured density of states divided by the BCS density of states above 11 
meV was not used in the fitting procedure that produced a^Fiuj) and a comparison of theory 
and experiment in the multiple-phonon region is a valid test of the theory. Reproduced from 
Ref. 1153 
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Figure 16: a'^F{iy) for Pb (solid curve) vs. u, along with the estimates obtained from Eq. ( p?] ) 
with an impurity scattering rate, 1/r = 1 meV (dotted) and 10 meV (dot-dashed). These 
are both qualitatively quite accurate, before they become negative at higher frequencies. 
Also plotted is the result (dashed curve, indiscernible from the solid curve) obtained from a 
full numerical inversion, as described in the text. Taken from the second reference in Ref. 
|T78 . 
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Figure 17: The a^F(z/) for K3C60 (solid curve) extracted from the reflectance data of Degiorgi 
et al. ||314|| , using Eq. (^Tf ). For purposes of analysis we have omitted the negative parts. The 
neutron scattering results from Ref. [|189|| (dashed curve) are also shown. Clearly the energy 
scale in a^F{v) matches that of the phonons, and some of the peaks even line up correctly. 
Finally, the dotted curve comes from an analysis of photoemission data [ p.9Cm , where we 
have arbitrarily broadened the phonon spectrum with Lorentzian lineshapes. Taken from 
the second reference in Ref. [|178| . 
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Figure 18: The universal function G{Q) as a function of normalized phonon energy Q = Q/T^ 
which enters the curve for the functional derivative of with respect to a'^F{uj) in the A®® 
model of Ref. 



195|, from which this figure was taken. 
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Figure 19: The constant c(/i*) in the relation /c^Tc = c{^*)A for the maximum associated 
with a given A as a function of ^* . Placed on the same figure (sohd dots) are the results 
for Tf./A obtained in the case of many strong coupling superconductors for which a'^F{(jj) is 
known from tunneling spectroscopy. The solid points all fall below the maximum curve as 



they must. Adapted from Ref. |196 . 
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Figure 20: The isotope coefficient, /3, vs. Tc, for various values of A and loe- Along each 
curve Tc changes because the Coulomb pseudopotential //* is being varied. These results 
show that a low value of f3 is difficult to attain with high Tc. On the other hand, for low Tc 
materials, it is not so difficult. 
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Figure 21: (a) The temperature dependence of the BCS order parameter, and (b) the result- 
ing densities of states at various temperatures below Tc. The only effect of finite temperatures 
on these latter curves is a reduced gap. 
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Figure 22: (a) A(ia;„) and Z{iu!n) vs a;„, the fermion Matsubara frequency, for various 
temperatures, as indicated. Note that the curves are relatively smooth and featureless, 
and at low temperatures little change occurs, except that more Matsubara frequencies are 
present. In (a) the units of A are meV. These were produced for Pb. 
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Figure 23: The (a) real and (b) imaginary parts of the gap function (in meV) on the real 
frequency axis, for Pb, for various temperatures, as in the previous figure. Note the consid- 
erable structure present on the real axis. Also shown is the (c) real and (d) imaginary part 
of the renormalization function, Z{u!) vs cu. 
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Figure 24: Calculated densities of states of Pb for various temperatures. In contrast to the 
BCS case (Fig. (21b), at high temperatures there is considerable smearing. 
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Figure 25: The ratio 2Aq/ ksTc vs Tc/uin- The sohd dots represent results from the full 
numerical solutions of the Eliashberg equations. Experiment tends to agree to within 10%. 
In increasing order of Tc/uJen, the dots correspond to the following systems: Ai, V, Ta, 
Sn, Ti, Tio,gBio,i, In, Nb (Butler), Nb (Arnold), ^35^(1), VsSi (Kihl.), Nb (Rowell), Mo, 
PboATio.6, La, V,Ga, NhAi{2), NhGe{2), Pbo.eTioA, Pb, NhAi{3), Pbo.sTio.2, Hg, Nb.Sn, 
PbogBio i, Nb-^Ai^l), NbsGe^l), P60.8-Bio.25 Pbo.rBio^, and P6o.65-Bio.35- The drawn curve 
corresponds to 2Ao/A;bTc = 3.53[1 + 12.5(Tc/u;in)'^^'n{uJen/2Tc)]. The insert shows results for 
different scaled a'^F{uj) spectra. They all correspond to the same value of Tc and of uen as 
Pb. They serve to show that some deviation from the general trend is possible. Reproduced 
from Ref. 
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Figure 26: Specific heat of aluminium as a function of temperature in the superconducting 
state and the normal state (applied field of 300 Gauss). Data taken from Ref. | p35[ . The 
BCS prediction, given the normal state data, is given by the solid curve. 
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Figure 27: The specific heat ratio, AC (Tc) / {'yTc) vs Tc/ugn- The dots represent results from 
the full numerical solutions of the Eliashberg equations. Experiment tends to agree to within 
10%. In increasing order of Tc/ujin, the dots correspond to the following systems: A£, V, 
Ta, Sn, Ti, T^.gSio.i, In, Nb (Butler), Nb (Arnold), V^Si 1, V^Si (Kihl.), Nb (Rowell), 
Mo, PboATio.e, La, V,Ga, iV63A£(2), NhGei2), Pbo^TioA, Pb, NhAii3), Pbo.sTio.2, Hg, 
NbsSn, Pbo.gBio,^, NbsAiil), Nb^Geil), P6o.85io.2, P&o.t^^o.s, and Pbo.65Bio.35- The drawn 
curve corresponds to AC(Tc)/'yTc = 1.43(1 + 53(Tc/co'^„)^£n(u;fe/3Tc)). Adapted from Ref. 

EH- 
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Figure 28: (a) cri(z/) vs. p in the zero temperature BCS superconducting state for the 
various impurity scattering rates indicated. The absorption onset at 2A(0) remains sharp 
independent of the scattering rate. A delta-function contribution (not shown) is also present 
at the origin, (b) Same as in (a) except for the frequency times the imaginary part of the 
conductivity. The optical gap is a little less evident in the dirty limit. The conductivity is 
given in units of ne^/m = Up/Air). Taken from Ref. |p.8H . 
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Figure 29: Frequency dependence of cri(z/) near the clean limit (1/t = 1 meV) for various 
temperatures in the BCS superconducting state. The appearance of a gap is evident, even 
at temperatures close to Tc. (b) Same as in (a), but for z/(T2(z/). The appearance of a gap is 
evident in the imaginary part of the conductivity as well. The conductivity is given in units 
of ne^/m = Up/An). Taken from Ref. ||181|| . 
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Figure 30: Frequency dependence of o"i(z/) near the dirty limit (l/r = 25 meV) for various 
temperatures in the BCS superconducting state. The appearance of a gap is evident, even 
at temperatures close to Tc. (b) Same as in (a), but for 1/(72(1^). The appearance of a gap is 
evident in the imaginary part of the conductivity as well. The conductivity is given in units 
of ne^/m = ujI^/A.'k). 
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Figure 31: Nuclear spin relaxation rate vs. reduced temperature. Data points for Indium 
are indicated by circles and triangles, while data for YBa2Cu307 are indicated by squares 
and crosses. The solid curves are calculated with Eliashberg theory for Indium (upper curve) 
and two model spectra with A = 1.66 and 3.2 (lowest curve). Agreement is good in the case 
of Indium and the lowest curve. Reproduced from [ p.68|| . 
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Figure 32: (a) Conductivity ratio, (Tis/(Ti„, versus reduced temperature, T/Tc, in the BCS 
limit, for various impurity scattering rates. From top to bottom the curves are calculated for 
= 100, 2, 1, 0.5, 0.1, 0.05, and 0. The frequency used was u/ Aq = 0.02. In the clean limit 
(l/r = 0) the coherence peak has disappeared, (b) Same quantity as in (a), but for different 
coupling strengths, A = 0.77, 1.5, 2.3, and 3.1. (The peak diminishes with increasing coupling 
strength). These were computed in the dirty limit (l/r = 500 meV) and for u — 0.05 meV. 
The result for A = 0.77 (largest maximum) is nearly identical with the BCS result, (c) 
Conductivity ratio versus frequency normalized to the zero temperature gap edge, u/Aq, for 
the same coupling strengths as in (b). The curves decrease in magnitude with increasing 
coupling strength. The maximum apparent in (b) for A = 0.77 and 1.5 is also clear here 
since the two uppermost curves have magnitude greater than unity. As the coupling strength 
increases the conductivity ratio becomes independent of frequency. Calculations are in the 
dirty limit, with T/Tc = 0.85. 
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Figure 33: Microwave conductivity normalized to the normal state, cris/(Ji„, as a function 
of reduced temperature T/Tg. The open squares are the data for Nb. The dotted curve 
is the BCS result with experimental frequency cu — 17 GHz and impurity scattering rate 
l/r = 100.0 meV (dirty limit). The solid and dashed curves are the results of full Eliash- 
berg calculations with two different {a'^F{u})) spectra. None of the theoretical curves can 
reproduce the data. 
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Figure 34: The real part (a,b) and the imaginary part (c,d) of the conductivity at essentially 
zero temperature (T/T^ = 0.3) with 1/r = 2 meV (a,c) and 1/r = 25 meV (b,d). In all 
cases we have used the BKBO spectrum scaled to give the designated value of, A, while Tc is 
held fixed at 29 K by adjusting ji*. Increased coupling strength suppresses both 0"i(z/) and 
va2{y) and broadens the minimum in the latter at 2 A. Note that 2 A increases slightly as 
the coupling strength is increased. The conductivity is given in units of ne^/m = uj^/Att). 
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Figure 35: Measured z/o"2(z/) vs. frequency at T = 9 K and at T = 300 K (solid curves). Also 
shown are the theoretical fits, using the BKBO spectrum, scaled so that A = 0.2 (dashed 
curves). is kept fixed to the experimental value with a negative /i*. Finally, theoretical 
fits are also shown with A = 1 (dotted curves). The latter curves are clearly incompatible 
with the experimental results. Adapted from Ref. ||186 
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Figure 36: The conductivity- derived scattering rate, 1/t(z/) = -^Re (l/cr(z/)) vs. frequency 
in the normal state for pure elastic scattering (dashed line), combined elastic and inelastic 
scattering (BKBO spectrum with A = 1), and pure inelastic scattering using a model spin 
fluctuation spectrum appropriate to YBCO. Because of the difference in spectral function 
frequency scales, the result for YBCO continues to rise with frequency, even at 300 meV. 



Reproduced from ||181|| . 
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Figure 37: Conductivity- derived scattering rate, 1/t{u) vs. frequency in the BCS s-wave 
superconducting state for (a) 1/r = 2 meV and (b) 1/r = 25 meV. An abrupt onset of 
absorption at the optical gap at temperatures near Tc is more apparent in (a) than in (b). 
The horizontal dashed line indicates the normal state result. Reproduced from ||181|| . 
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Figure 38: The conductivity-derived scattering rate, 1/t(z/) vs. frequency in the s-wave 
superconducting state for (a) 1/r — meV and (b) l/r = 25 meV, for temperatures as 
indicated. In both cases we used the BKBO spectrum with A = 1. In (a) there is no 
signature for a gap, while one remains at low temperatures in (b). 
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Figure 39: The widths of low energy [(^(^0]Ti acoustic phonons broaden appreciably at tem- 
peratures above Tc, the superconducting transition temperature. This figure shows the same 
phonon profile above and below Tc ~ 18.0 K Reproduced from [305 . 
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Figure 40: (a) Real and (b) Imaginary part of AIl{u + iS)/N{0) vs z//(2Ao) at zero temper- 
ature, for various impurity scattering rates, l/(rAo) = (solid), 1 (dotted), and 6 (dashed), 
in the weak coupling (BCS) approximation. Below twice the gap edge the phonons soften; 
above twice the gap edge they harden in the clean limit and soften in the dirty limit. Note 
the narrowing that occurs below the gap edge in the presence of impurity scattering. 
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